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1. DOUBLE-LEVEL METALLIZATION TECHNIQUES 

Most of the work during the first quarter was directed toward the 
development of a sputtering system for preparing aluminum and aluminum-alloy 
films. A photograph of the completed system Is shown In Figure 1.1. Briefly, 
the system consists of the following. 

The basic vacuum system Is a Varlan-NRC 6 Inch oil diffusion and mech- 
anical roughing pump equipped for automatic and manual operation. The 
diffusion pump Is equipped for LN 2 cooling of the cold trap. The sputter gun 
and power supply were obtained from Sloan Technology. The sputtering chamber 
was designed and built at Mississippi State. It consists of a Corning 12X18 
glass cylinder and an aluminum top plate machined to accomodate up to three 
sputter guns. A cold cathode discharge gauge was constructed and installed 
in the baseplate of the chamber to measure the pressure during the sputtering 
operation. A throttle valve with several threaded holes for accomodating plugs 
Is operated by one of two mechanical feedthroughs in the baseplate. A lift 
mechanism with a reversible motor was designed and constructed for raising 
the top plate and sputter guns. 

inside the chamber Is equipped with a rotating table which accomodates 
up to eight wafers of 1i-2” In diameter. The table Is driven by a vacuum sealed 
shaded pole motor through a magnetic coupling at 7 rpm. The entire motor- 
table assembly Is rotated by a chaln-sprocket-mechanlcal feedthrough arrangement 
through three sputter-gun and two mask positions. The mask Is attached to the 
rotating assembly and provides one hole through which the sputter-gun deposits 
metal on the wafers. A crystal film thicknes sensor Is located beneath the 
sputter-gun and receives a deposit through a hole at an unused wafer position 
on the table. 



Figure 1.1 Photograph of metallization system 
using sputter-gun source. 
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The sputtering chamber opens Into a class 100 clean bench In order to 
maintain a high level of cleanliness. The system Is located In the metallizing 
and bonding room of the microelectronics laboratories In SImrall Engineering 
Building, and this room was designed with air conditioning and filtering units 
to maintain a class 30,000 environment. The entire system has been constructed 
and checked out and Is ready for depositing sputtered films. 

The Installation of a six-tube Thermco Ranger diffusion furnace was 
completed with the addition of a venting system for exhausting the scavenger 
boxes. All that remains to be done Is to line the tube and c»nnect the nitrogen 
ambient source In order to anneal the aluminum films. 

2. TW0-Dlf€NSI0NAL MODELS FOR MOS TRANSISTORS 

The work done during the first quarter was directed toward preliminary 
investigations of numerical schemes and computational algorithms for solving 
the semi conduction equations for a two-dimensional field. 

A recent report has described the application of the finite-element method 
to the analysis of a JFET.' The finite-element method has been used for some 
time in solving problems In mechanics and elasticity; however, it has only 
recently been applied to semi conduct I on problems. This method has the power 
to treat some problems, such as eigen-value problems, for which the finite- 
difference method Is awkward If at all applicable. It can also be applied to 
the solution of field distributions governed by partial differential equations, 
and one of the most attractive features as compared to the finite-element method 
Is purported to be the esse of treating non-rectangular geometries and Irregular 
boundaries. For example, the geometry of the VMOS structure could be accom- 
odated. It was decided to further investigate this technique. 
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In order to better understend the eppi Icebf f Ity of the method, It was 
appITed to a one- dimensional linear diffusion problem. This simple problem 
Is one for which familiar results are available for comparison and at the 
same time taxes the finite-element method. In Its most valid form, the finite 
element method is applicable to variational problems in which a true minimum 
of an energy-related function exists. Such a minimum does not apply for the 
semiconductor problem In which current flow occurs by diffusive and conductive 
mechanisms, it has been proposed that a "weak form", the so-called Galerkin 
method, be applied to such problems.^ The typical semiconductor problem Is 
a non-linear boundary value and initial condition problem of which the linear 
diffusion problem is a very special case. In the example chosen, the diffusion 
variable, u, obeys: 


U ( X ^ t } 1 

* u^ » 1 

(a) 

(2.1) 

U^(x,t)|^.^ 

• 0 

(b) 



■ 1 , x»0 
0 , x>0 

(c) 


c 

• 

c 

X 

X 

II 

f(x,t) « 0 

0<x<a 

(2.2) 


The Galerkin formulation of this problem is: 

,a 

/q ( UfV - u^Vj^ - fv ) dx ■ 0, (2.3) 

where v(x,t) represents a "triai function" which Is used to approximate u(x,t). 

The finite element methods uses a set of "hill functions" as illustrated In 
Figure 2.1 to construct the v(x,t) approx I met Ion. Two of the popular hill 
functions are the Hermlte bicubic and the bilinear functions which are Illustrated 
in the figure and were used in the example. The final form of the approximate 
solution is; 
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N 

v(x,t) • ^E^qj(t) ♦j(x) . (2.4) 

On the node points the solution Is approximated by the set (qi(t)) for 
the type of hill functions which overlap as illustrated In Tlgure 2.1. 

Solution for the set (q|(t)} Is then analogous to solving for the set (u|(t)) 
on the node points using the finite difference technique. The equations for 
the set (q|(t)) are obtained by substituting (2.4) Into (2.5). 

N a 

£ /o ♦(♦j ♦ <i|-iiii4 ♦ ♦♦!> • 0 

l«i at •' ax ax'* 

J « 1,2,3, N (2.5) 

From (2.5) a set of time differential equations Is obtained which is 
solved using an implicit numerical method. 

The solution of the problem posed by the example is closely appro! ximated 
by the erfe function In the range 0<,x£3 Ifa»6, and this solution 
was used to compere the accuracy of the finite element and finite difference 
methods. Figures 2. 2-2. 4 shew the maximum error as a function of the recip- 
rocal of the number of grid points and the size of the time step, l.e. 

At/Ax^. 


The error obtained In the solution by the bilinear finite element 
method Is very nearly the sme as that obtained with the finite difference 
method. This was not surprising because the systwn of equations for {q|(t)} 
and (u|(t)) were quite similar. What was surprising was that the Hermite 
bicubic finite element produced such poor results, although this surprize 
was based upon the Intuition that since this element was more difficult to 
uso it should provide some reward for the difficulty. 
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This experience resulted In some skepticism that the finite element method 
would be effective for semiconductor problems. A more recent paper has rein- 
forced this attitude.^ This paper Indicates that the proper formulations of 
the semiconductor problem for the finite element approach remain to be demon- 
strated, and. In agreement with our observations, point out thAt the application 
of Galerkin’s method Is subject to skepticism. Therefore, It was concluded 
that the further work should be based upon the finite-difference method 
which we have used before although the finite-element method Is Intriguing 
and may be further developed in the future. 

The second phase of this program Is underway to develop a computational 
program with which to generate data from a model. During the past quarter, 
the rtajor emphasis was upon deriving a simple two-dimensional algorithm which 
could be economically used In a simulation program during Its development. 

It Is plausible that this algorithm can be later modified to become adequate 
for modeling of other effects which are significant In certain situations, 
e.g., short channel effects, avalanche breakdown, etc. 

During the past quarter an algorithm has been developed based upon the 
usual assumption that the mobile carriers are Included In an Infinitesimally 
thick layer of charge at the Sl-Si02 Interface. The current flow is then 
described by a one-dimensional equation: 

I * kT u.,( + quiE-S ) , 0<x<a (2.6) 

r ^ ^ ^ 

where I and S are channel current and charge per unit channel width, Is 

the tangential Interface field and and are the mobility factors 

accounting respectively for gate modulation and hot electron effects. This 
equation Is solved Iteratively with Poisson's equation which includes two- 
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dimensional effects: 

■ -p/e. (2.7) 

Equation (2.7) will be solved for one segment of a periodic structure with 
respect to x and equation (2.6) will be integrated to produce auxiliary 
equations for the boundary conditions at the Si-Si02 interface. It will 
be assumed Initially that S ® 0 at x * a, the point at which the normal 
component of the interface field changes sign. The solution algorithm then 
proceeds in an iterative fashion to solve (2.6) and (2.7) simultaneously 
for the potential distribution from which the current, i, is ultimately 
calculated as a function of the gate, drain and body voltage . 

The model at this point admittedly has some short-comings, mainly due to 
the neglectlon of generation-recombination mechanisms. Therefore, It will 
not treat Impact avalanche and bulk generated leakage currents. It Is believed 
at this time that such effects can be treated by adding another iterative loop 
to the algorithm. The major emphasis at this point will be the development 
of a program for input and output data management and including subroutines 
which generate data internally within the program and solve systems of equations 
which will be encountered In implementing the algorithm. 
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3. REDISTRIBUTION DIFFUSIONS FOR I ON- IMPLANTED PREDEPOSITS OF BORON AND 
PHOSPHORUS IN SOS FILMS. 


The objective of this work was to produce curves describing the variation 
with diffusion time and temperature of the Junction depth, sheet resistance 
and Integrated Impurity dose. This data has been generated for boron and 
phosphorus redistributed in nitrogen, dry oxygon and steam ambients for <11 I> 
oriented SOS films. The following section presents discussions of the Implant- 
ation and redistribution model, further program develop, the computational 
procedure and of the ccxnputed results. 

3.1 The Redistribution Model: 


There are three aspects of the redistribution model which are considered: 

(a) the Implanted protllo, (b) the oxidation model, and (c) the diffusivlty 
model . 

(a) The Implanted profile. 

The Implanted profile Is described by the gausslan function, 

y , (3.1) 

L Rp 

where C Is the concentration, y Is the distance from the entrant silicon surface, 
Rp Is the range and AR^ Is the straggle for the Ifi^lant. The peak ccncentra- 
tlon, Is related to the implant dose, Oj^p by: 


Redistribution data has been generated for the following conditions: 


^Imp* 

5x1o’^ , 10’^ 

, 5x10^^ , 1o'^ cm"^ 

Rp: 

0.2735 Min 

80 keV boron Implant. 

0.1727 " 

150 keV phosphorus. 

ARp: 

0.0665 pm 

80 keV boron Implant. 

0.0440 " 

150 keV phosphorus. 
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ORIGINAL k’AGLjfc 
OF POOR OUALITY 

The doses are light to moderate resulting In concentrations no heavier than 
18 —3 

6x10 cm , and the range-straggle values are typical of those employed at MSFC. 
It Is assumed that all of the Ions become activated shortly after redistrib- 
ution begins and thereby diffuse by a substitutional mechanism Involving 
vacancies. 


(b) Oxidation model: 

The oxidation model is assumed to be the same as for bulk silicon and 

5 

the data of Deal et. al. has been used to calculate the oxidation rate accord- 
ing to: 


dx^ = B/ ( 2xo B/C ), (3.3) 

dt 

where B and C follow Boltzmann- 1 ike temperature dependences. Figures (3.1) 
and (3.2) illustrate the oxide thickness dependence upon time and temperature 
for both dry O 2 and steam ambients. 

During the oxidation, the silicon film thickness is reduced according 
to: 

W * Wq - oXq, (3.4) 

where is the Initial film thickness, taken to be 1 pm, and a = 0.45 

is the ratio of the densities of SIO 2 to silicon. Redistribution data Is 

o 

given for = 1 pm and an initial oxide thickness of Xq = 300 A. 

(c) Dlffuslvlty model : 

The dlffuslvlty model for boron was discussed In an earlier report® and 
it includes a linear dependence of the diffusivity upon the vacancy concen- 
tration as well as the field-enhancement effect. The dlffuslvlty model for 
phosphorus Includes only the field-enhancement effect which is sufficient to 
describe the non-linear behavior of phosphorus diffusions at concentrations 
lower than lo'^ cm”^ as shown by Barry^ and Fair and Tsai®. The dlffuslvlty- 


ORIGIN RANE is 

OF POOR QUALITY 
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Q a 

temperature dependence is after Fair and Fair and Tsai adaptation of data 
by Ghostagore'®. For either boron or phosphorous the effective diffusivlty 
Is given by: 

= D(u) X { 1 + u/ / u^ + 1 ) , (3.4) 

where, 

u « C/ 2 nj, (3.5) 

and, 

D(u) a Dg u, for boron, (3.6) 

a , for phosphorus. 

and where nj is the Intrinsic carrier concentration at the diffusion temperature 
and D| and Of are the Intrinsic diffusivities of boron and phosphorus: 

Og a 3.17 exp ( -3.59eV/ kgT) cm^/sec. , 

D« a 3.85 exp (-3.66eV/ kgT ) (3.7) 

3.2 Further Program Development: 

The program which was used to generate the data has been described 
In detail in an earlier report. It was noted that the program was developed 
in such a way that one could take advantage of a normalization procedure for 
predeposition diffusions and generate data applicable to different film 
thicknesses. However, It is not possible to gain such an advantage for re- 
distribution diffusions Involving Ion-implants or growth of an oxide. Then 
the program was used to generate data, it was discovered that some other features 
of the program are extraneous unless further refined. 

The program was developed to account for both thin and thick oxides such 
as would be encountered in some practical situations. However, such a simu- 
lation requires the incorporation of a warped grid system, a modification which 
would require considerably more effort. Therefore, the variable oxide feature 
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Is of limited value at this tlrne, since the program, at best, only approximates 
the conditions for growth of a very thin oxide during redistribution. 

A modification was made which allows accurate treatment of redistribution 
under oxidizing conditions when only a single oxide thickness Is Involved. 

The original program treated the oxidation process with regard to the boundary 
conditions; however, unlike the case of bulk silicon, one must also account 
for the reduction of the silicon film thickness. This feature Is now included 
In the program. During the simulation of a redistribution in an oxidizing 
ambient, the vertical grid spacing continuously shrinks while the horizontal 
grid spacing is constant. The modification does not show up on logic flow 
diagrams at the level of detail which has previously been given. For complete- 
ness, a new listing of the affected main and sub-programs is given in the 
appendix. 

3.3 Computational Procedure; 

The program described in an earlier report, and modified as described In 
the preceding, was used to generate the data. Two-dimensional data was obtained 
In the form of Isoconcentration contours for typical situations. The bulk 
of the data which can be correlated with experimental measurements Is generated 
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using a quasi -one- dimensional model in a manner described in a previous report. 

A brief review of the procedure is given in the following. 

For generation of sheet resistance, junction depth and Integrated Impurity 
dose data as a function of time and temperature, only a one- dimensional 
profile need be calculated. This Is accomplished by making the horizontal grid 
only three units wide but keeping the field six film thicknesses wide. 

Periodic boundary conditions for the horizontal dimensions are employed In 
the program and result in a calculation which produces the vertical profile 
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equivalent to a none-dlmenslonl model. Thus without sacrificing the generality 
of the program for treating two-dimensional cases, the amount of computing 
time Is drastically reduced when the data that Is desired docs not require the 
full power of the program. 

The vertical grid varies from thirty one to sixty one points as required 
for accuracy In details of the profile, and most of the data is not sensitive 
to the number of grid points used If the number is chosen In this range. For 
the purpose of Illustrating the unusual nature of phosphorous profiles, the 
larger number of points was required. 

3.4 Discussion of Results: 

First, some of the unusual behavior of redistribution diffusions In SOS 
films will be discussed in this section. Next, tie format for the calculated 
curves will be discussed, and, finally, the bulk of the generated data Is given 
In the appendix without further comment. 

Figures (3.1) and (3.2) Illustrate the oxide thickness growth and silicon 

film thickness reduction as functions of time for <111> silicon films oxidized 

In steam and dry 0. ambients. The evolution traced ^ginning with an initial 
o ^ 

oxide of 300 A thickness on an SOS film of 1 pm Initial thickness. The curves 
are shwon for four temperatures. The data are necessary for interpreting some 
of the results for simulated redistributions. 

Figures (3.3) and (3.4) show impurity profiles for boron and phosphorus 
implants being redistributed in a steam ambient at 1000 deg. C. The profiles 
are all plotted with a common origin as would be the cose for experimentally 
derived profiles where the SI-SIO 2 Interface would serve as the logical origin. 
However, the profiles are normalized with respect to the film thickness which 
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is of course shrinking. The boron profiles are not unusual but show the well- 
known leaching effect due to segregation Into the oxide. The phosphorus 
profiles show the effect of Ifi^urltles being rejected from the oxide. There 
is a pile-up of Impurities In front of the advancing SI-SIO^ Interface and 
then a dip which eventually disappears. It Is easy for one to draw an erron- 
eous conclusion from observing the profiles, because It appears that the 
Integrated dose should increase for at least remain constant and the shoot 
resistance should decrease with tlfr«. This Is not true. Although the segre- 
gation coefficient favors phosphorus in silicon vs. 3102* eventually all of 
the phosphorus will be in the SI02 when the SOS film Is cumpletety oxidized 
since the model assumes that there is no diffusion Into the sapphire. 

Figures (3. 5-5. 7) illustrate the behavior of the junction migration, 
sheet resistance variation, and Integrated impurity dose variation rver a 
long period of time. All of the curves are plotted with respect to normalized 
time, and true time Is obtained by multiplying by the normal izeng time value 
given on the plot. Junction depths are In microns, sheet resistance values 
are In oivns, and dose values are In cm~^ units unless otherwise marked. The 
curves are given In the typical former for all of the data. 

For an Ion-Implanted profile, there are In fact two Junctions until one 
of the junctions emerges at the SI-SIO2 Interface. Therefore, the sheet resist- 
ance values are for the burfod layer until the front Junction disappears. 

This typically occurs In a short time compared with that for through-diffusion 
of the back Junction. Figure (3.5) Illustrates the through-diffusion of the 
back Junction for the hMvIer doses. This always occurs for redistribution In 
a nitrogen ambient but not necessarily so for an oxidizing ambient. After the 
through diffusion, or even before for light doses, the Junction depth will 
eventually decrease due to the reduction of the film thickness or due to th? 
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relatively slow advancement of the Junction with respect to the movinq 
SI-SIO 2 Interface. In some of the data presented In the appendix, the junction 
appears to remain almost stationary for this same reason. The variation of 
the sheet resistance and dose with redistribution time also may appear strange 
when compared with results for bulk silicon; however, consideration of the 
previously mentioned factors also explains these results. 

Two-dimensional Isoconcentration contours are given In the appendix for 
the various ambients and the two Impurity types. The results are not as 
remarkable as those given In the last report which were for chemically pre- 
deposited boron. In that case, there was initially a heavy concentration of 
fast diffusing Impurities at the Sl-Si02 Interface which were strongly retarded 
due to the segregation phencxnena. This does not happen with the lon-lmplanted 
predeposit because the Initial profile lies below the interface at which the 
segregation phenomena Is eff-. 'tlve. 


26 


REFERENCES 


1. J.J. Barnes t A Two-D?mens tonal Simulation of MESFETS . Technical Report, 

RAX-TR-76-r53, May 1976. 

2. G. Strang and G.J. Fix, An Analysis of the Finite Element *1ethod. 

Prentice-Hall, 1973. 

3. J.H. Hohl, "Variational Principles for Semiconductor Device Modeling 

with Finite Elements," IBM J. Res. Develop., Vol. 22, Ho. 2, 

Mar. 1978, p159. 

4. D. Wdo, Private communication provided the range-straggle data. 

5. B.E. Deal and A.S. Grove, "General Relationship for the Thermal 

Oxidation of Silicon," J. Appl. Physics, Vol. 36, No. 12, 

Dec. 1965, p 3770, Deal, Sklar, Grove and Snow, "Characteristics 
of the Surface-State Charge of Thermally Oxidized Silicon," 

J. Electrochem. Soc., Vol. 114, No. 3, Mar. 1967, p 266. 

6. J.D. Gassaway and 0* Mahmood, Computer Simulation of Two-Dimensional 

Impurity Diffusion In SHI con-on-Sapph I re FI Ims . Interim report, 

NASA Contract NAS8-26749, March 1973. 

7. M.L. Barry, "Doped Oxides as Diffusion Sources," J. Electrochem. Soc., 

Vol. 117, No. 11, Nov. 1970, p1405. 

■ i; • , ■ 

8. R.B. Fair and J.C.C. Tsai, "A Quantitative Model for the Diffusion 

of Phosphorus In Silicon and the Emitter Dip Effect," J. of Electrochem. 
Soc., Vol. 124, No. 7, July 1977, p1107. 

9. R.B. Fair, "Boron Diffusion In Silicon-Concentration and Orientation 

Dependence, Background Effects, and Profile Estimation," J. Electrochem. 
Soc., Vol, 122, No. 6, June 1975, p 800. 


APPENDIX 


BORON DATA 



o o) CO r- Ln 

OJ 'a t ^4 4 >. • I « 4 

lU t>l iu LU LlJ UJ 
O ». ) O CD O O 

« -4 % I ... I 

I I i I I I 

□ O < + X o 

o 

00 . 

: -J' CD CD 

CD CD CD 

■ CD • 

CD - CD CD 
II II II II 

LU 

rc: 

D n.. 

I *- 1 i.j 
CL I • 

rc (Ji 

LU 

a_ i u lu 
>r; >:: >: 

01 lU ‘ ^ 

X I - L •• L - 


P/VUL 



C 3 — ^\J <n sf 

■ • • ■ ■ 

o o O CD CD 



original 
OF POGR 



3.3 




I 

I 


r 


E; 

f 

I 




o tn CO to 

OJ •> -( w 4 *-4 mi m-i 

Ul LM Ui LU Ui LU 

o c:> o o o cj 


I I I I I I 

H O ^ 4 - X ^ 
CD 

03 . 

: C D CD 

CD CD (DJ 
. CD CD . 
CD - •■• CD 
II II II II 

I J 

rc 
3 n.. 

f - I lJ 
( i: I • 

rx: CO 

hJ 

a. hJ I d 
>: >D 

oj I L.l * 


NITROGEN AMBIENT 










0E20 


ORIGINAL page J© 
OF POOR QUALITY 


tn 05 u) 

%-4 VI «4 **l 

III lU UJ III UJ 

<5 O O O O 


I I I I I I 

El O < f X O 
CD 

CO . 

; :P CD CD 
CD CD 
• CD CD . 
CD V--I CNJ CD 


Id S 

nr; 

. j n.. 

I - lu 

CL I • g 

cc 0 5 5 

(i.j S 

a. I iJ Ld * 

>D ic g 

CM 1. ij ► -H P — I K 











^ 









t 









/ 

f 









t 










































— 

0-- 











1 




-X ' 


-X 





\ 






1 



I 






i 



1 






i 

1 









i 

i 

i 









\ 

4 

C 


O 






1 









1 

5 

i 

I 









1 

1 

j 







0.0 0.2 0.« 0.6 0,8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 3.4 3.6 3.8 


UJ <fJ L<J UJ Ul UJ 

O Cl c* o o o 




I I I I I I 

a o < + X 
o 

oo . 

: CD CD 

CD CD CO 
- CD CD . 
CD OO CD 
II II II II 


LlJ 

rc 

:d 0- 

I - lU 


cc cn 

iij 

n_ i.ij LiJ 

>D >1 s: 

CM liJ • 


ITROGEN AMBIENT 







ORlCmAL PAOW 

OF POOR QUAUTx 


A 5 





t 


/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 


f 


O CO 00 LD 

t\| I V « 

tu i>j Ui ui III uj 

O Cl o o o o 


-m 4 *-4 m -4 » 4 »<4 

I I I t I I 

H O 4 + X O 
CD 

CD . 

: :i^ CD CD 
C D C D CD 
. CD C D . 
CD '-H in • 
II II II II 


a 


lU 

re: 

. D o„ 

I - I u 
a; I - 
nc CO 
lU 

re I u ijj 
>D >I >Z 
oj I-lJ ^ 

X I -• h -I— 


t; 


L 


Nitrogen ahbient 
















NCRHALIZING TIME >6716.0 HIN 
fl»S:EN‘T « NITRCGEN BORON 


ORIGINAL PAGE IS 
OF POOR QUALITY 


nin* n r-r 


‘nrrnT'nrn**"^ 


irnrrr I n' 


€Nt 


i I 



Si 




C) 


o 


o 


4 . ’ 
f ' 


cj 

rc> 




CJ 


C > 

Li* 


O 




v.j 

Oi 


J J — I .llliU .UL^ ULUJ-LLJ _ J 


-03ft)r^iD4n <«> <\i 


0 

1 


r9ion*0 in ^ m 


ct 


^•<ll3r><Oi© ;i* ir> iv 


o 


I 

o 


KJ 

CJ 


o 


UJ 
>: 
► 4 

H 

CJ 

UJ 

M 

fc> 4 

. J 

u: 

cr: 

t.) 


33NUlSrS3y J33HS 




A 9 


0 
CD 

CO tu 

1 « 


o 


«.> 


UJ 


o i; ♦ I 

o I: 

o 

z 

«-» » 
hNl 

►- 

« -I 5 f. 
CC Mi 

0. x: 

z rn 

tti 13 ». 
9 - Z Cfi 


iin rr r i — rrmTO'*! — oi i n n ~ 


1 

m 

*o 

s 

“t 






O 

S 


(W 


/ 


.LIJLlXi-J I UILIX U .J, 


-U.UJJ,..L,L^ 


GO 
> •< 
o 


j-OA^r-Oin sf m 


nj 


fv 

o 


*o>car'^cdtr> z o» M 


CO 

o 


Z ID cv 


c:» 


*. ♦ 


c. • 
r- 


o 


«■ 

<n 


o 


c • 
ir> 


o 


o 

: J* 


KO 






• ) 
V t 


O 




l/> 

» * 

o 


Ui 
>■ 
^ 4 

I ‘ 

n 

UJ 

IVJ 

I* I 

. I 

d: 

>: 

cr 

*“^d» 


JUiyndHi Q3iyy93iNi 


NORMAL 



Hidja NOiiaNnr 






jLiiyndwi a3iuyo3iNi 





3.0Nt)lSIS:^y J3.1HS 


TEK? = 1000. C 

NORMRLIZING TIME =<4114. 1 MIN 

RX3! ENT = NITROGEN BORON 


ORIGINAL ?AGR IS 

OF POOR QUALITY 



A 15 



xiiundwi adiuyodiNi 


EMP = 1050.0 THICKSE 


original page is 

OP POOR QUALITY 



HJd3a NOIiONfir 















0E20 




cr* CO in 

• » I V •£ » < •* 4 

UJ i 4 J UJ Lrr lij 

C) O O O O 


I..LJ 

cc 
:) 0 ^ 
i ~ ijj 

a: h- 

(x; <n I 

a. ijj ii.1 ^ 

cj I.jlJ • •-■■Hfc 

< ^ L I _ = 




1 h' 



CM 

CO 


lo 

mrntm m9tm^ 

iO 1 

• 

• 

• 

• 

• 

• 

o 

o 

o 

o 

o 

O ( 


« 










O CO t“- <n if, 

C\J • . » %--K • 4 «. I « 4 

UJ i* l«J IgU iU lu 

O < 1 O Cl I. . O 

• *4 r 4 ,4 , 4 

t 1 t t i I 

o o ^ -I* X o 


CD 

cT3 

CD 

«'D . 

i 

C D CD 

CD 

CD O 

C\J 

• CD 

CD cn 

C D -4 

CO : f 

il il 

II II 


IlJ 

(T. 

-■) n_ 

f ” lU 

a: \ - 
cc: 01 s 

UJ I 

a. lu UJ 

^ 

UJ * 

I - , i g 





- -K v. . 



E2C 




E20 


A 

•> 


ORIGINAL PAGE IS 
OF POOR QUALITY 

A 24 








TEMP * 900 



:3DNbiSIS3y 13DHS 


N3RMRLIZE0 TIME 


ORIGINAL PAGE Ib 
OF POOR QUALITY 



Jii]yndw[ a3iduo3iNi 


NGRMqLIZED TIME 






TEMP = 950.0 

NORMALIZING TIME »483.3 MIN 

AMBIENT = DRY OXYGEN BORON 








'rrrrrT 


TTr 



L 111 .L-LJ — L UJJ.U...I.J L 




OJ 


n 


-ilLUJ.. 


CO w 


»-cwDr-c£>4/^ ^ <n 


01 


in 

o 

^4 


CO 

CJ 


fVJ 

o 


33NHiSIS3H J.33HS 



A 30 


LU 

5: 

^ i 

h- 

iL\ 

hsj 


cc 

5: 

rr 


TE/.P = 1000. C THICKNESS = 0.0 Ci<- 

NORMnLIZING TIME -120.0 HIN 

RKeiENT = DRY OXYGE:'i BORON 


ORIGINAL PAGE IS 
OF POOR QUAUIY 



HldBQ NOIiJNnr 


CdznbwyDN 



TEHP = 1000.0 

NORWRLIZING TIME =120.0 MIN 



TEMP = ICOC.O 

NORMaLIZING TIKE =120.0 MIN 

RMBIENT = OP.T SXYGEN BOB0N 


ORIGINAL PAGL i- 
OF POOR QUALITY 



Ai]yndw[ a3itiyo:jiNi 




TEMP « 1050.0 

NCRMRLIZING TIME -30.0 KIN 

RMSIENT - CRY CXTGEK BOROM 


ORIGINAL PAGE IS 
OF POOR QUALITY 


rmTT "I — r 


"rnrnrr 


Tn- 


'ft|i 


n“TTnr“T — i 


A 35 


O 

O 



IJXLLU-JL 

cs 


to 

o 


J. 


— HU-LLX.1 

o 

10^(0 

O 


« r^»lrf J-J [mm L 

M ^otr^ioxn M 

o o 


o 

ID 

o 


:a3Nyisis3y x33hs 


.SO C.^iC 0.50 O.oO 0.70 C.80 C.&O 

NCRKP.LIZED TIME 





0*0 






ORIGINAL PAGE IS 
OF POOR QUALITY 



o 




f 



zr 

CO 

CO 

CO 

ft 

<n 

rJ^ 

O') 

<M 

ft 

o 

CO 

00 

ft 

CVJ 

a> 

ft 

OJ 

C\J 

OJ 

ro 

o 

« 

Od 

CO 

cD 

=J- 

rj 

o 

00 

d 

<D 

ft 

c» 

d 

Od 

o 

o 

d 


A 39 


o 





TEMP » 900.0 THICKNESS 

NORMALIZING TIME « 1400.0 MIN 



Hid3Q NDIlONnr 


NGR^RLIZEG "riMF 







NORMALIZING TIME 



NORMflirZED TIME 






TEMP * 950.0 

NORMfiLIZING TIME =483.3 MIN 



:33NbiSIS3y 133HS 


.to C.2C 0.30 0.40 0.50 0 

NSRyiRLIZED TIKE 












TEMP » 1000.0 

NORMRuIZING TIKE =120.0 MIN 
fiHBIEM = ETEPK BORON 


ORIGINAL PAGE Ifci 
OF POOR QUAUTY 



iLLLLLU LUILU X J LU IU.U,J, 


=J» 

o 


*WS <OLn:t> CO CM 


CO 

o 


fO rsi 


ry 

«-4 

o 


.•wvana* fo cm 


"irriTT'rrT' 


.Jil.LLLU_i 

9 

CO CM 

•4 

o 


juiyndwi aaiduoaiNi 










TEHP » 1050.0 

NORMALIZING TIME -30.0 


r 



uiyndwi oaxyyodiNi 


JGRMRLIZE- 


APPENDIX 


PHOSPHORUS DATA 




















> 


B 4 



t 

O 


(\J CD 

• • i 

o o o 


ID to r- 

• • • 

o o o 


QD 

O 


o> 

o 


I 




w 


u 


okigwal 

OP p(.K)H QUAlilTY 


B 5 



L 




O Q) CO n- to LO 
«>« 

Ul UJ LU UJ UJ LU 
O O O O O CD 


I I i I I I 

0 O <1 + X O 

oo 

o . 

ooo o 
oo o 

.00 . 

O— *LO^- 
II II II II 

LU 

cc 

O CL_ 
f— LU 
CL l~ 

QI CD 
LU 

Q_ LU UJ 

s: s: s: 

<y LU * — H-H 
X h- h- h- 



2.8 3.0 3.2 



PAGE 



I 

o 


Hld3Q NOriONni 



TE«P = 900.0 

NORMALIZING TIME »1UC0.0 HIN 

R«3I£«T = NITRCGEN PHOSPHOROUS 



jmyndwi oaibyoaxNi 


NGRM3LIZED TIM 




950.0 







THICKNES 



grmoliz 


















TEMP = 1050.0 

NORMALIZING TIME «30.0 MIN 

AMBIENT = NITROGEN PHOSPHOROUS 


original 
np POOR QUAUi' 


OF POOR 


B 18 


inTn~i‘'T'’T 


'rm'i'TT'-T n mn 


t 


o 

•V. 

II 


.IJ 

% 


9 


H 

<n 


W 

% 


<n 


o 

CJ 


o 
— J cn 




1LLLLXJL.J LUXLLLX.-J. 


-IIUJJ^L. 


in 

o 

•»>4 


ron^r^oinr* (o 


<\j 


a* <n 

S’ 




CO 

♦-M 

CD 


r<Mr^(nKO :r <n 


01 


OJ 


o 

Cl* 


o 


o 

r- 


o 


o 

CD 


o 


o 

lO 


o 


c> 


o 

OJ 


CD 


o 

ru 


o 


o 


o 


o 

o 


CD 


XJliyndWI Q31HU031NI 


CR^^flLIZED TIME 






































TEMP . 900. C 

NORMALIZING TIME -mOO.O MIN 

RM3IENT = DRY OXYGEN’ PHOSPHOROUS 



;33NtllS[S3U 133HS 


NORKRLIZEO 



TEMP « 900.0 

NORMftLIZING TIME *11100. C HIN 

RKSIEN7 * CRY OXTGc*.' PHOSPI 




nn* 


L®, 




■■rr"T — inTTTT'T n 




T 


ITTTTT r"l ■ 


5 




H-UJ-LJ L iinui 1 1 


ilUJ.X,.! .1 L 


to 


zr 

o 




« sJ* *n A4 

<o 


fM 

*4 

o 


,-CKr-<OI/> to 


Ci 

c .1 


o 


cr» 


• 


o 


o 


CJ 


» 




O 


r- 


• 


C5 




fi:? 


rj 

UJ 

>: 

C-» 

V -1 

» - 

UT/ 


* 


O 

rj 

UJ 


M 

o 


rP 







ir. 


o 


if 

« 


O 


c.» 


C**i 


• 


c* 


*.■> 


• • 


o 


o 


o 


• 


<r» 



O 

^4 


XUUOdHI a3iyH03.lNl 



IUd3Q NOILONfir 




TEMP = 950.0 

NORMALIZING TIKE =183.3 MIN 

AMBIENT = DRY OXYGEN PHOSPHOROUS 




mTn"n — 


“IIT 


TrrT'^Ti 


TTTTTTf 



LiLLLLU — L LUJJLU-.J I L- 1 , J I 5 

Cl o o • 

CO c\j ^0m'^coio=j' m w -acrcou) co cm o 

in a< CO ry 

C) O o C) 

«-4 


; 33 Nt)lSI 9 ^y 1 -J 3 HB 


N3RMRLIZED 




uiyndN[ a3iyy93XNi 


NQRMRLEZED TEHE 


THICKNESS 



ORMflL 







TEMP = 1000 



JlilHndWt Q3ibH03iNl 


SnyiPLIZEC? 






ORHRLIZI 








TEMP = 1050.0 

NORMALIZING TIME =30.0 MIN 

fiMSIENT = DRY CXY6EN PHOSPHOROUS 



XllHOdWt a3i«9031NI 


ORMflLIZED TIKE 











CE20 














o <x> r< 

*-M » < ^ 9^ * 4 9*4 

ut i.«i Ui uj tit t«i 

o o o o o 

• •*•»• 

1 » < 

I I I I i I 

PI 0 < f X <> 
CD 

O C3 

O . 

O O O 

• O CD (DO 
• O CD CO 
O ■> * < ' CO 
II (i II II 

lU 

It: 

D n„ 

i - 1 1 J 

a: » - 
nc: CO 
ui 

UL. 1.U l.U 

>:: >z >z 

«M I lJ ' -< • •< 
X I - I - I - 











QtfAunj 



HJd3n NOiiONnr 






TEM! 

NOR! 



jLiiyndwi a3idy93iNi 





ORIGINAL PAGE it 
OF POOR QUALITY 



o o 


Hid3Q NDIiaNfir 
















TEKF » 1000,0 

NORKflLIZING TIME *120.0 MIN 

RMdIENT -- STER.^ PHOSPHOROUS 


OKIUINALPAO^ 
OF POOR QUALmt 


B 



X ^I1±±.JLJ Ill- UJ 1 I I 

:* <D ru ^mo'^win:^ <n m sf tn ot 

) o O 


o 

o 

o 


Jii]ynciw[ Q3ibiyo3iNi 


.20 0.3G O.IIC 0.5C 0.6C 0.70 O.SC C.SO I . CC 

N3R.«PLIZED TIME 





B $3 



O 

o 


o 


o 

Cw 


Cj 

I' 


o 


CO 


o 


r 3 - 


o 


o 


o 

CVJ 


Ct 


o 


o 

o 


C) 


o 


CD 


c> 


uiyndwf a3iUH03iNi 


\’0R?<r5'.IZE3 TIME 


^fC» 0 N#SnS 2 il IfSHPi oT 

C ••• SHPirJN ^InT F»tOr,HAli 


I 

7 

7 

•I 

I 

iS 

11 

12 

\l 

IS 

tf 

II 

I* 

11 
21 
in 

I! 

if 

21 

2f 

!? 

U 

32 

3<l 

it 

37 

it 

10 

11 

12 

13 

10 

14 
^7 

4 i 

If 

10 

11 

4 S 

It 

4 1 
Sf 

I? 

42 

43 

44 

45 
44 
47 
41 
4f 

5 ? 

H 

74 

74 

77 

74 

74 


••• D 4?4 orrv ••• 

Clio 

SW .-.8 «, 

r ICLO II 10 _ 

• C«IC| 4 TM ANft *Th C«I 0 

jcoD > luT n roi <uO svlo 


Tint 2 l 0 ti 


jcoto 

14210 

NTff4 

NKfTn 


Or olv. 
or olv. 


IN 

In 


YNI 4 V 
YW 47 V 
iniNV 
1 N 4 XV 
ritio 

0l«CN4I0N 


SUtrsNlD OfV* Y 
1 OtYtS |A SuoOO. 
fUOftllD 0|Vt fN % 
i OlYfS I SuilOD 
niH Y4k IN 
NM VAt fi^ 

NtN VAt IN 
NA« YAk IN 


4X{« 
4«U 
AXIS 


AXIS 


n I 0 t 4 cid«s 

OsllOt 


vu 

*>IS 

«l|s 

(KIS 


. iroi , iJi 10 , I ool ,oi 10 , luQi «T |HC< io» ISO* 

oiMr»<oiON iiiiOtioO) ■ 

OfMCN^ION |K(S>,tvTS|,|LI2l,XO(|n,IOOt 
••• ICAO IN 0«T« 

• r *0 lOOfN^i T,NST«*I ,N^T^ 2 ,„ 8 Ta 3 ,i,i$T^«I,NSTP$ 
so 10 l•l,N|»LT 
irii.ro.n hst^»», 8 Tpi 
ir(i«ro» 2 ) 

irii*co<ii N$T^»N 8 l ^2 

iril.KltNI uST^-NfTP'. 

IflI.rS.Sl MSTPa*iSTP 8 

no 10 J>I,N 8 TP 

Nr«OM**tlONl jN«itiU*>lNT 

'•l 2 h I{'?^.t"A*| 9 CLT,Drtt.r 01 tT 

I, loll •I* f .Ji , 1 ] u I ji , x J i t . ji ,r 

II U »•IOlT C' 
n <2|*'V6 Cn 


IC 


•irAOMN 
ICAOM" „ 

iri i«mont.cq,i ) 
ir< i«MiNr.co,i I 


|0< I . Jl •TI"CI I • Jl ,101 I • J) 


|F( l*NINT,{g,2l 
IM UN0NT,t(i,2l 
ir < i«MiNT,ro, 3 i 
in iahont,Co,3* 


f) ; I l••N|TllOC• 
ll t2)a*CN • 


I0« rotHArioiiOi 

101 rONMATI |H0«4C|S**I 

102 rO«M«TI2(8**l ,7| in* 
IC2 rONMATINI lO.NflO.A* 

10** rOtNATI iHOtOI lOl 
I0)> rPNH«TI2A*) 

<•09 rOllMATI |H 0 < I Ol • * XxlN , 

• //lHO,3i,<Mrio,S.'«<<>i 


,s* .•»«», 


, 5 * ‘TNIN • »,Sl *TM*I 


OATi m 6T ,lN«x,yHt|/0.!i*7&,7. ,8, / 
THAXaTMil/An, 

ICalO 

• •• IN|T| aTC THt »*l oT 

no II Nat, 3 

»£Ao 102 . iini*,i;i,A*,«iYm.iai.ii.jcno 

• CAO lt'» Iftlm, Jr.'*l0,**7Sfl»,Nl5Ta,ffl|ny,»H 
FAINT <IOO,XM|NV,|aAlv,TN|Nv,TNAl« 
iriJCnO.Nl.ni tcnna2 
irijcoo.lOtOi lcnn*2 
fAU. AioTSt I0,0,in«9| 

CALL ALOr 1 |•S«|.C•-3l 

• •• RAAA NOAOf A 

.'ALL ALOT|0.0,TM«(,2..I 1 
TALL AL0TlAMAi,VMAl,2.*l I 
CALL flOT |lMAI.0.n.2,.n 
CALL FtOT (3, 0,0,0. 2,-1 * 

OtA* 6AI0 

V 0 lVa»N»I/ri niTl.lAHlDi 
J«a I 

00 12 |a|,Jr.H|0 
nn 13 J« I . I( ,NTSTn 
VSFCaALOOlOl A lOAT I J* I I afn I V 


VHaBv.XMinV.INAIW 




•0 

• I 

:! 

II 

tt 


CALL AlOT.O.O.TS^rtAl 


li 


CALC ^lOT(C*l2iV<tfCt2.ei 
CALL ^L0T|lMAItTt»C»|.0| 


C2 


o*Nisr» 


>Si»iiuI/Ki«tit>»ieii 

VAttSvMlltV 


irooi 


ORiGriML Pack 

OF Pooh qi’autv 


SB Ti 

bO 17 J*l»ie»NT 
ir(itCBtj«»ei*A 
vtPC«AbO«iOift.e 

Ut51H!tf .. .. 

Htrietiiii! 


i*ANO.J*Nt*U «e TO 

ft.eAT(J*j«i»*roiv 

aoo 


IT 


•0 fo 17 ‘ ••OA*VAt,«..»l I 

aoo ♦o»io.io.'.o»i»i i 

|7 


VAL i•VAtl•lf<« 

I* JR«JR«iO 

IValMAIV/rbOATI IftRlOt 
l6tl01«l6Rln«| 

liTCat i-l 

itTei«i$Te*fi»2 
VAL.aXV*M*ll 
i* CALL NUMRCBiISTCI 

• •• 


fit L.tlLi - - • 

?.‘tt s?W8t!i;s8::5;i;;i;ijini;;8:;54, 
jsti !!S8jk!|:j:}:i:s{:;{SJ:,:, !!•-’' 


Sf*t 
MUM 

CRLi. SVMI 
• 0*0. )2I 

rALl MUMI 
CALL «f«i _. 
CAU. STMtOLt 
iriN.ce 




•21 CALL RrH»Ol lR*l«k*i*M«T** 


• 0* • 1 0) 


• 0 B 221 -1»**«T,»TmIC«NC$$ • 

trfA.re*2i cAtu MiiMBrtis.A,2.k,Mi;r.7oisT,o*ii i 
••• ORA* CORVtS 

TlOCaVRAl/riOATt jeRlOi 
IV|alH«l/AM«AV 
VHWRai ,/TH|MV 
50 22 I»I.RRLT 
(ALi RLOT(0,0«0*5.2I 
iMi.ri.it nstr>n*!m 
IMI.C 0.2I NSTR«NR7f2 
lMI.eO.2l NtTR«MRlR2 
IMI.IO.RI hSTR«MiTR«| 
lMI.rO.it tiSTRaniRTRi 
no 21 ja|,h«TR 

iru.ro.it t.a) 
tru.Ne.it uaj 

lltOVC* • I NC I I # J I • I “ I 
|MN.e0.2t 60 TO 7 CI 
tr(N.ro.2t 60 TO 7P| 

VHOVtaTL06*AI.06l^|RS| i.jtavNVRt 
, 65 TO 2t 

20 1 VMeveafL0b«tk06IB|0l I . Jt aTMVN t 

. 60 TO 21 

202 VNOVean06*tLOCIO|R.M I . JtavMVRt 

28 (ALL RLOT(laOVt.rMO«r,L.Ot 


8 »In», 


I 


iriN.Me. 

(ALL RLO 
50 22 Ja|,N«TR 


>2t.6 

iTtO 




40 

iru.co.li L>3 

C5 

41 

42 

irij.NCtii l>2 
XMOVCbT|NC( 1 • 


43 

4<l 

XI ( 1 • Jl>l{ ( 1 • j)*lvO|ST~0«*l»»XO( 1 .JM/I 


4S 

44 

47 

4t 

44 

VMOVCbVL06*AL061O|X| II |J1*YHVR) 
23 CALL PLOT|XmOVE»VnOVc,L«6) 

22 eONTINUr 

CALL RL0T(O*iO*««**> 

II CONTINUE 


70 

71 

STOP 

FNO 



C4 


l|C*0N*SnS2l I I aPARtM 

1 subroutine P«RAHlrl«T| 

i IHRLICIT double NbECIsION fA-H»0-2l 

* 2 fall j|{*7*ie*IO**«RTlH 

B r6»C6/lBt*2e*S*(TI ) 

10 r!«!i?B$el4*MT)*«l t*Bl 


u 

n 


«ANo'.R*fiT»l *OoS) so TO 2 

return 


ORIGnMT. PAGE IS 
OF P(X)R QUALI'ni 




MCR 0 N*SnS 2 f 1 I tOXDATA 


C6 


1 

2 

% 

il 

II 

11 

I? 

22 

23 
2 H 
2S 

24 

12 

Zf 


5U8»0tjTlNf OX0 *T*,AMbnT,o»INT,T,».c,H,KB| 


SUnffOUTlNC ^08 OX|nAT|ON ^ARtMCTeRS ANO SCGKCCaTION CfflfHt 

*•••••••••••••*•••••••••••••••••••••••••••••••••••••••••••••••••« 

IM^UCIT OOUILC ^»CCMl0NlA-H,0-n 
iNTeeK* ahbnTiOpmnt . 

OOUBie PRCCISION KB«N,N|tMa 

iruHHNT.eB. I » eo TO ii 
B><t**t 0277 o-IO«Dex»<~*«^S* 78 /ri 
IMOBlNTtCQ. I ) Ca4*V<iBSBD 


12 




triOBlNT.CQ. n Ca4*V<ISSBD-i«0CXPI-22lB«l.o7/Tl 
60 TO H 

ssriimTo •yaeCXP 1*1 3«lA(Ait<tf /Ti 
iriOB|NT>EQ. I ) Ca2i00f J0-|a0CXPI>2‘i| IBtVn/r) 
ir(0BtNT»CQ.2) Can t 332770-1 aorXPf.BISSl •*8/T) 
IF I OR I NT >£0.3 I Ca7*09aRSO-|aOEXF|-2RBS7.o2B/Ti 
CONTINUE 

N|a33.3*0EXP(«0.87/<KA*T| I 
M3a20,O*0EXP(>0«<%7/(Ka*T) I 
irtORINTtEO.I I NaMi 
IFIOR|NTiEO< 2I Na(Nl«M2)/2,0 
IF(0R|NT«e0»3l HaH3 

return 

END 


ORIGINAL PAGE IS 
OP POOR QUALITY 


^ICKONaSnSZdurNONT 

SUlitOUTlNt rSONTl , LM.it ,TtMP»OrI (OtLT.JsTtP, 

• «type,koti,*ohtJ 
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so 


*0 


sua. Foa calculatimG contour front novencnT data and 
NCA0IN6 AND HR|T|nr Thf SaNE ON FILE 

implicit OOUSLE PRCClRlONlA-H.O-n 
OIhEnSION XFRoNTiaI* Y^NonTU) 
data lOOl.lBLNK.tRTARdHt.lH ,|M«/ 

IF(KTVPc«NE<I I SO TO rOo 

Mead file 

DO AO KK;I|JSTEP 
READ ITiSTOI mark 
IF iNAMKtNEtlOOLl On TO SO 
R|A 0 J»iN 30 l TFMP.nPI.oELT 
“»3T0» “ 


RCAOlTi 


mark 


READ (T.R30) (OuM.LsI.AI 
READ (T.Sanl (OiiM.Laj.A) 
REAOIfiNZO) TmAXI ,«(Ma|I .K.lM 


CONTINUE 
, RETURN 
aoo CONTINUE 

00 2A0 LL-l.A , ^ 

CONVAL> 10 . 0 **FLOATl 70 -LL«n 

XFbCONDEPI IHAXI ,JMaxIiOi*JNaX| .O iO»CONvALi 
TFaCONDEPI IMAXI . JmAX | . <>2 ,6 .0 .0 tCOMVAL ) 

IF iXFtEO.O.Ol «0 TO.2S0 

XFR0NTILLI«( AF-fLOat(KI l•IXolST/FLOATI t MAX | -I I I 
2S0 CONTINUE 

IF (TF*EO.O.O) eO Tn ZAO 

VFRONTILL t■IFLOtT UmAXI i.TFl/FI OAT UMAX 1-1 I 
2A0 CONTINUE 

STORE contour front MOVEMENT DaTa |F IFILf > I 
HARKalOOL 

IF (KOTI.gT.I) HARK-IrLNK 
WRITE 19.3*01 Mark 
IF (KOTl.GT.ll GO. To 2^0 
RRITCI9.**30( TEHR.OFI .OELT 


MARKbIBLNK 
RR|TCI9.390i n«KK 
270 rritE (9.930) (XFrONT(LL) .lL"I tA) 


WRITE (9.9301 <YFrONT(LLIiLL>I.AI 
WR|TE(9.920) |Max I t JM aAI .K.LN 
390 FORMAT ( IHQ. AA) 

900 format (In ,3fI0.9> 

920 format ( IM ,9|5) 

930 format (IH .A(C|9.9.2x) I 

return 

END 
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46An IN r«OH nATA OCCtC 
jsrcp* 4 or TINC pTCPP to oc pcao 


IN 


ITCHP • SINUUatIOn tCMPtPATUnr*# 

•on - I anoa*«2«i 

•OLT • tINC 5UP«* 

OOUOLC PPCCISION tir40liTf400i4l •Tr40NT(40f1i4} 
T€HP«0rifOLT 


400I 




200 


I 


OIMCNSION 
eiNCNSION 
OATA lOOLt 

o*ta ia/Ihc _ 

OATA IX/AHIN X O.AHiRrCTI tAHON 
... DATA tt/AHlN t O.AHlRrCri tAHON 
390 rOKMATHHOiAAl 
*•20 FORMAT! IH i**lfc» . . 

••30 FORMAT! IH iA!Ct*!.9i2xt > 

••• RCAO DATA FROM nECX 
RCAO 200i JSTeF 
FORMaT!OIIO» 

T|MC«0* 

00 I XK"liJSTcP 
ReAO!B«39&T mark 

usstS'iASirffW.iif.isJ 

RCAO!8t390) HaRK . 

REA 0 !Bi*( 30 i ! XFRQNT!KK«tLI .LL«I .41 
READ I A. *• 30 ) ! VFR 0 i|T!KK«LU.!.L>i .41 
REAO!0.*!20» IMAXI .JMaiI.LM 
TINCaTIHC^Dl T 
TAUlKKt>T|Me 

continue . 

DATA H&TiXMAX.TNAI^0.n>7S,7*0.S.n/ 
••• initiate THE Pint 
00 00 N-1.2 . . . 

CALL PLOTS! 10.0. I4*0l 
CALL PLOT! I.S.I.K,-3i 
••• ORAM SOROfR 

CALL PL 0 T! 0 . 0 .TMaIi 2 .» 1 I 
call PL 0 TIXMAXiTMAXi 2 .>|| 

CALL PL0TIXMAXi0.ni2««t> 

CALL PLOT(O.O.O.Q, 2 .>i I 

••• scale grid 

TTaU»TAU! jSTEPl*inO» 

IO|V»INT!TTaU) 

IS»IOI V/2S*25*25 . 

TAUMaFLOATI iSl/IH'i* 

VTIaXPRONT! 1 . 4 ) 

IP!N.E0«2) TT|>YFR0NT( 1 .41 
00 2 I«2.JSTEP. 

VT2>XPR0NT ! I .41 
|F!N.C0>2) YT2*YF40NTI I .4) 
YYY«AMAXI!YYI .YY?I 
2 YYI"YYY 

VAP«YYY*tCO. 

IY0»INT|YAP| , 

IT»IY0/25*25»25 
YOMAXbFLOATi iTI/irO. 

••• ORAN GRIO 
OXaXHAX/20« 

00 10 I*! I 19 

«»FLOaT ! I I •OX 

CALL PL 0 T(X, 0 .t 3 l . 

CALL PL0T!X,YMAX,7,H 
lO CONTINUE 
OYaYHAX/20. 

00 20 Ja | • 19 
YaFLOAT IJI ‘OY 
CALL PLOT ( 0.0 ,Y , 3 ) 

CALL PLOT ( XmAX i Y , ? • I » 

20 CONTINUE 
a»a erase LABIE ARFA 
00 30 Ia 3 ,l*t 

XXa ( I • I ) aOX 

call PLOTIXI, | 7 .anY, 3 i 
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*7 
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|00 

181 

lOS 

|04 

|07 

loa 

|07 

l!I 

1 u 

lU 

\\i 

lit 


c 


c 


AO CALL ALOT(XX.|7*&«OV,|2| 

00 AO |■lf•70 

ty«M-i »*of 

CALL ALOTIO«.VV,31 
AO CALL ALOTt lA.tDl.yYi |»» 

••• AOT *YH|Ol 

CALL SYH|0Ll0.K»A.7,e.>AI>*tA»0.e«S3l 
iriN.rOtii i«»iY , . 

CALL $YH80L< l•A•«t.AS,n.t^|^. ix.O.OilAl 
••• AXES NUHSrAS 

88ft8Al? Iil 

^ VALBriOATlt.||«TAuM/|0* 

&0 CALL NUNaeR(X(Y,HfiT.VALt0«0»2> 

X ••0 • s 

00 AO Jali2| 

VbYLOATI J»n»0Y-,«4J 
VAL«rL0AT|J«| |*YOMAX/»Ot 
40 CALL NUMBtAlX,Y.H«T,WAitiO»Ot22, . ^ , 

call SYM60L(-«7fi,>< i.iAISiToiSTanCE' f«0* »•! 

CALL SYMeOLC2.»-.%t«lYl3t*»!OAMALIZtO T|Hr • TAU»,n. 
00 70 LL«I»4_ _ 

CALL AL0YI0.0»0*0»3> 

IM0VC«0« , . 

CALL AL0T(TAUtXH0VCf2.0l 
00 70 ML*1»JSTEA 


TAUA»TAUI«L|*XHAX Jtaum^ 
XMOVC*XFRONYIML«LL> •vHAX/YOMaX 


IF (N.E0*2» vN0VE«yER0MT(NL tLL t*YMAX/VOMAX 
IF<TAU.6 T.O.O*ANO*VHOvC*EQ* 0*0*AMB«N,ES*7I iO TO 70 
tFIN*E0*2| XHOVEaVAOvr^ 

|FCXhOVE*lT.O> I XMOVCaOtC 
CALL ALOTlTAUFf XMftVE,7iO> 

70 CONTINUE 

. CALL ALOTlO. »0» .Y4*» 

00 CONTINUE 
ATOF 
FNO 


• 21 > 
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10 
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n 


suBdouTiNi SMriLr •TiHr*0£tT»oti-T.^eH^,TM»K,M( 
.JH*m .TOISTiXO.UMuNTI 


Sub WKtres secr BrsiSTtNci, 

iNPURirr ON UNIT M 


junction OCPTm ANo 


INPLICIT OOUBLl PbCCIrIONI t 
NRITCI 13*2001 JHAXliUPRNT 

***** 58|Ttcli!I?6ol '?e«P'*TH*K,0ej.T,OCLV,T0IST 

_ «BITCM3*I00) RS«TJI *vJUNC.O*TlMr *X0 
too rORNATt iHOiAClS.ti 
RETURN 
END 


TJI .TjUNC, 

intiorateo 
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i|CI»0N#Sft5?l I ) .OUTPUT 
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7 
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1 1 


IS 

II 

il 

I? 


2S 

2 * 

27 

2 i 

2 f 

I? 

li 
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s* 

37 

31 


C 

C 

c 


SUKKOIITiNE OUTPUTlXiV, |HA«I| JHAk I 
•IO«IT|HCiXO,PTtHE,DT|MCt lAHlNTI 

t.ansiint data p.inT ouT 


KiLH»JjiT|MriT 0 |ITi 


INAllitllHAlU 


implicit DOUILC PtEC|«lONlA*H.O»7l 
COMMON /CON/CIt ( AfiiMl 
OIHCNSION KIDiTlitilntl.liaOIII 
IMAXalMAXI-l 
|MAX3a|MAXUI 

pSiSl loo! il!l|Mc|PT|MC.OT|MCiK,XlK).l 

PMINT lOti INiNat. iMAXl lil 
print 102i (Xh I .t>2,tMAXI,2) 

PRINT lOS 

8«T0|ST-T»JI 

IF( I AMRNTiNP» 3) aBUM«X|.Jt*4l*/CL0ATcJMtX|>| )) 
2 PRINT I03i Oi ICBI I 1 1 JA; l*2, IMAXI .21 
PRINT lOti lX0(nila2.|MAX|t2| 
irMAM0NT,Nr.3> PfiNT jOR# «| 
io» pormati/iho.iox,*: si.riLM • • *.riOtSi 
*|0X|*« C“ ‘ 


lOR format* iho 


OX|nC THICKNCSR IN CM • * / 1 HOi 1 31 • II ( 1 Pr 10* 3 1 1 


PRINT ION, JJ 

too rORMAT<lHO>IOXil2wTtHr,STCP > • |N«3X,7NTtMC • iFin«3. 

•SX,*ClAPSeO T|Me ,N SFC<t*,2X,<PRtDCP > • •FI0*3|2II I *OR|Vr IN • «, 
• F|0«3// 

•inx.ioxioc POSITION*/ 

•lox.txi • ii2. • i • .leA.x.ax.txit.ii,*) • ••c*t2//i 

tot FORMATI/ZilHO. )HI • .AX.|2|I0| 

102 FORMAT ( / . I HO > 3 X • $HX a, i 7 X . I 2 < I PF 1 0* 3 1 ) 

103 rORNATdH, 2X,*V«itlPr2.1,7X,llnFCI0.3|| 

ION rORNAT(//,lHO,'NO. Or ITERATION • ••ISi 
I05 FORMAT! IHO) 

10* FOMMAT! |H| , lOXt l5A‘*iTR«,*T|MC»tAAI 

return 

FNO 


PRT S0SF.PL0T»C0NT0UR,.SUBI0N<.TR|nA(i..AnCt.XTZ,.PL0T..COM0CP.«MAlH 
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.PUOT»CONTOg« 

SUatOuTiNt PLTCON»a,CONVaU,NC,U.j«» 
KOCONCCNTKaT ION #iOT aUlPNOCOaN r 


Cl? 


8IKSII8S ««:8::<?!}???4?:j?!}?8?:8!4W1??!1. 

DATA 2EN0/l.0C*2nA . .. 

RATA l(CSl|,JitJ«lj2t,l>|.2l/atl«-l*0,0.«.»0«l|0*OiO«»/ . 

2 * 1 * t il^n 1 * I •** 8 • »C* 

0*0**i *0 tOjOi *1 t 0«0*p« { *0* I *0 *Otn I* I >0 lOtOr , 


fi*I* ! or /o*p*o*Oii*o.i*Bto* 

OATA l(0SI|.jlij«ti2t,l>|*t)/O*»*l*0 - - - 
DATA IMaIi VNAI/0.a*N.a/ 


i« 6 * 6 «d**( 6 * 


• 0 * 0 / 


c* 


NCMIsNC*! 

eONTURaALOClOICON\|AL| 

• NlttT PEN T0,0»»«|N 

SCAkXaXHAX*l •k/| |v*2) 
8CAkY«»NA*/( JE«I I 

• STaNT CONTOUN 9PARCM 


10 


c* 


c» 


NT >0 

IL«tK-l 

00 fO 

•‘^locate i^UANE caOSSiNCf 

•yTyU*2k»Uk06IO(AI I i JII^ALOGlOl Al l«l * Jl I «Al06l 01 A t | 1 1* 

•ALOSIOUI 1*1 I J*l I 1 1 
all lalO^OavAl I I 
aiai«ai I l 

• li'CaTe TaiANCura 

00 0 Ka| |N 

mPal 

60 70 121 .27.23, ?4tl .K 
ai2l>Al 1*1 .Jl 
aiStaAl I ,JI 
60 TO 30 
ai2>aAl I .Jl 
alSIaAl I t J*l I 
60 TO 30 
al2)*Al I .J*l I 
aniaAl 1*1 .J*| I 
60 TO 30 
ai2l>Al 1*1 .J*t I 
ai3laA( 1*1 .Jl 
60 TO 30 

• locate INTENSECTlONa 

00 38 Hal, 3 

ir ICONVAL.LT. AN|mMR(H) ,B(M*| l l . 0 a.C 0 NV«L* 6 T.AMAXl IBIHI ,a(M*| I I I 
• 60 TO 38 

«iPaNP*l 

aaaAL 06 IOIBlH*l } | .ALOaiOia I M | I 

ir (Aasiaai .GT.zraOl «o to 33 

60°!^ 3** 

TXaOSII.HI*cS(l,HI*D 

TTaO5l2,H)*rSl2,Ml*0 

«INT*| .NPlaoTI l,»»*CT( I ,l,al*TA*fT( I . 2 ,Kt*TV*|l 
T|MT*I .NP|aOT(2,K |*CT(2.l,«l*TX*rTI2.2.K|*TT^JJ 
CONTINUE 

IP INP.LE. I ) 60 TO NO 
NTaNT*l 
NVINTlawP 
CONTINUL 
CONTINUE 

• scale points 

IP (NT.ce.Ol 60 Tn 8C 
00 A8 Kal.NT 
NMaNVtKt 
00 AS L«I.Nm _ 
k(ir«L}«x(KtU«SCAi X 
V (IC ,L I u *l ) •SCAI 7 
fOKTiNUE 

• PLOT CONTnUP 


20 


21 


22 


23 


2 ^ 


30 


3 S 


^0 

SO 


4C 


4b 



cn 


eo 71 K«|,Nt 

So 7|®ti|;A8******’’ TIK, II. OtlfliNCm 

CALL ^LOTIair.LI.vU.i I.2| 

n eoNTiSul***’ 

.. • HOVC ^CN TO OAftlN 
CND 
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<4tCtON*$nS2(ll*Trt|OA(i 
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li 

!l 

If 


fUBROUTiNt TB|0AB«JF«| I 


IMfLIClT OOUBLC RbCCIBION |«*M« 0 *rl 
rOMNON /TBI/ AURI tl(*f I ,C|«f I lOUf I .VUbI 

OIHCNSIOM tARMAIABltBrTAUf I 
Bf TBIJf jr» , , 

BBHHAlJf ^■D| jr l/UrtAj .if I 

jrBlajf ♦! 

DO 100 |•JrBl•L . 

tAtT«i-jr 
VILIaBAHHAUI 
00 200 K«IiLAST 

|al.oK 

too V( I lafiAHHAl I l-CI I l«V( t*l l/lirT«( I I 
BrtUBN 
BNO 
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•i|CR0N*SnSI( 1 1 • suit OK 

I SUBOOUTiNC .r.yOIST,rOtOP| 

1 §****Sgi5ooT*Sc*roS*«lSr«jT*N«*ION*tN^*J2T*olTj****** 

€••••••••••••••••••••••••••••••••••••#•••••••••••••#< 

IMPLICIT DOUILC 0»IC|«IOMU*M,0*n 
common /CON/CMI UlltBMI 
OlMCNStON VIII 


IfTvijMlJtT.I.O-Ji 


istfi 


.00 . 

• MONON IMMlANT ttSViOolo* RfV * oti t 'CST aM • ••CI|o*l 
DO I 0 JM*I ’ . . 

1 roi I I oJ>*C^***Of *K**OoO«l lrOIST»T|J|«MMI/OMMl»*ll 
00 2 J»l oJMtXI 
, DO > l•^oR. 

1 CDII I oJI«Cl| I I oJt 

2 rONVlNUC 
MCTUMN 
CNO 


• mICMOn'/ 

• |TM 


M|CM0N»SnS2l I loADC 
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J \ 


4 

7 
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f 

10 

11 

12 

\i 

li 


It 

I* 

20 

li 
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2 * 
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I! 

y% 

34 


fUAHOUTf N| Ate I tHilKl I J^AKI oX«V»CM«T|Hr fKoOri otCN^.rUAIoKnTt »XOo 
• rotsTof ahanti 


• • 

c«« 

c*»« 


SUttOUTINC MRIUS TMANtUNT DATA AN OaTa rlU ON UNIT i| 


VCO 


200 


I 

ICO 

JCO 


• • 
• • 
• • 


DOUtir MMrCtSION ftl oI»V,T|M* .TcmMoTNAIoAVI oXOiVOIST 

rOMNON /CON/CtI I A moANi 

DlMCNStON XOI I I 

DlMCNtlON I I OVII I 

data lOOLoltLNK/Utoloo / 

koti«coti*i 

NARKalOOL 

triKOTtoOTol I M«Rr*l«i MX 
* MAtX 


KMircti I oVOPi ) 
rOtMATl IHQ iAAI 
triKOTIoSTol I I 


80 |0 2 

MRITCIll ilOOi or I .ThxxoTcnAoToIST 
•mrcni o20puiinaxi ..imaxIoIanonv 

rOtMATI IN oti iOl 
«Rirci!l olOOl IXit oI>2oIM«X{| 
tiopi Ivnl olaioJNxxh 


NRITCIli 


NAMK«]IILNK 
•Rtfr i IT oKoni 


MAAv 




iMf 


rONTtNUC 
MRITCIll o200i 
NRITr I I I 0 tool 

S8,(r/TOSi!‘i;i,u 

RRlTClII olOOl ixniil 

rORMxTIlH otClSoRl 

MRINT 300i ROfi 

rORNxTIlHOtlOXo'KnTi 

RETURN 

rND 


,l» ol 

fit 


l> 2 o iMAxi I 
l-AXI I 


IIO/I 



Clfi 

<t|CllOM*ifiS2| I I.IV2 

I tUnNOUT IMC IT2I l««*l .jNAXI .«• T .in.TIMriK.Or I ,TtHp,TMAl,iori ,in, 

! c««*****I?!n*«»«!l*I«* ••• •••• •••• •#•••• ••• •• •• 

I C tuMOuTlMl RCAOt fuC THANSICKT 0«T« flLC nN Un|T h 

f * oeuliLt ^AfcisiON fti tll▼•T|Mc,rcMA•tN«»,nri •>e«Toisf 

■ rOMMOm /CON/r«| 

ftIHCdSION «o7l I 


MCAOn I t200t Kj I 

i88 

KC«0(|l tlOOl 

■“**“■**• •-»- I V < • • . I - ■ . IH* • I I 


roitr 

lAMtMT 


KrAoiiitiooi 
•CAOMl.lOOl 
CONT INUC 
ffC«OI 1 1 tlOOl 
RCAOiUtiooi TiNr 
• jajMAII (I *•! 


(Tm »l*t iJNAll I 

hark 


i 

KOO 


00 I J«JMAA| (I 
RCAOni tiOOl icti 1 1 .J| (laS, IHili ) 
R|«OMt|IO0t IIOI|>tt«2itH«Ril 

rORRRTi |HOi<A» 
iRRir RCRieor 




ARRlT ptaiovi''. 
|M«l2a|HAI|«l 
|H«I a|NAA|>| 

DO 7 ja|,jM«X| 

COM I *JI*COM)iJt ^ 
roll IMAAX.JiacIl ( t**AI.J) 

rcturn 

(NO 


.plot 

SURROuT INC PLOT (CRi lHt>l • JMAX I ,*u 


suRROuTiNc Plots tmC t«o dimcnsional propilc 
In thc outtput printout if iplot ■ I 


OOUBLF PNCCISION CBI.rS 
COHNON /CON/CBI URfSRt. 

OlPCNSION SVMbLIFI I fUi iNCIRSt 

(iHRllH tlHCilN iIHOiIm (IHEiIh iINF.IH • 

•|HG«|H •|nIi|M «|MK/ 

ORTA 00T,ST«R/lH,.tH*/tKP/2l/ 

OAT* |X/6«lH • IhViIm t lH«t | hX« IHI • IHS«tO«Ih / 

1HAXb|NAX|-I 
PRINT II . 

«T«ALOGIOl2.2nl t t 

pACTa|ALOGlOUSI>tTI/|FLOAT(KP>| I I 
K»0 

no 2 jajNAXI.Ifl 
00 I l•2.|MAXt 
IFICBI I I t J) .LC.O.n) GO TO 3 
, K-( (AL0GI0(CB| ( I , |l »- aTi/F*CTI*| ,0 
J CONTINUE 

IFIK.lT.II ULlNC(t>*00T 

IF (K*6E« It ANO.K.L p.KPI ULINCC I )•SVMBL (K I 

1 IFIK.GT.KPI ULINFiII-rTaR 

2 PR|'Nt''!o!7VlLI| (IiLINf ( I I •l■^t IMa^ i t lULtNC ( I ) I i•INAK.2t•l I 
PRINT 12 

PRINT 13 

PRINT INtSYMBLlKP) iCS 
KDaKP-l 

00 R I-XO.I.-I 
CBL»|Ot**(l l-l )»F4CTt.»T) 

CBHalOt**! UFACTtiTI 
R PRINT ISt SYMRLI I I tCnmCBL 
PRINT 16 

PROFh C»/ 

,»|M0.RI < IHUI ,///! , 

l| rORMATlIH tlX,50nH*»/<lH^,lSXi»X AXIS*! 

|3 fORNaTI /tIHOi’SvMBOLiC REPRESENTATION nF CONCENTRATION RANGfSa 
•//IHO.T 10. •SYHBOl.lt T3 *i» ’concentration RAnGC'/|H .TlO.61 Im>I .T30, 
•|RIIM->/ 'HO.TI3,|H*,T29,tAB0VE RURFACE rONC.tl 
|R FORHaHIH .TI3|AI «T29.* at fS • •.ElO.31 

FORMATIIH .TlS.At ,3X,.LCSS Th AN t , e | 0. 3 . 3X t • GRE *T eR ThAN nR EGUAi T 


PROFh E»/ 


|R FORHAT I IH 
FORMAT! IH 
^•0* .f 10*31 
I* FORMaTIIM 
RETURN 
ENO 


,Tl3.lH..T29t«BEl0R l.o30EII*l 
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.COnoEP 

rUNCTlON C0 n0c'*(M.N,1 .J,N|N,HaX .rf'NVALJ 
• locates concentration contours aLONO CiTMER 
A VCRTICAl OR horizontal 6H|0 ( iNE* EITHER 
LINEAR OR L06R|tH|N|C INVERSE INTERROLaT |0n 
CAN 6E USED. 


■EiNft CONTOURED IDINENSiONEd 
iN calling RROgRAN 


I.J 


CONVAL* 


• MIN, - 

• Ma* 

• 

• CONOCP- 


NON'ZERO VaL.UE SREClEtrS SitID line 
TO BE c0NToyR|0,_RoS|T|VE VALUE EOR 
linear interpolation. NECATIVr POR 
LOGRIThIOHIC interpolation. LiTHER 
I or j HUST^BE zero* 
cONOUR Value 

MINIMUM AND maximum S1IRSCR|PTr OP 

GRID l;NE tO be contoured. IHaV BE 
zero tr ENTIRE grid LINE IS TO BE 
contoured. I 

POSITION Or CONVAL On gR|0 LINE* IP 

conval is nUT OP Range of gRIo line 
VALUES conoEp Returns a value op 

ZERO. 
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implicit double PrEC|sI0N(a*H,0»ZI 
COMMON /C0N/A(AR.aRI 
no too I«I»N 

00 100 JMltN 

too |FIA(I(JIaLT*I* 0<«00| A( I ,jt>| ,0-100 
iLOGal 

IF l|*NE,0> GO To to 
IINCmI 
iMINaMIN 
|NAXaMAX-| 

IP ( MlN.EO.O) iMlNaj 
IF I MAX. EC. 01 iMAXaH.t 
IF IJ.LT.O) ILOGa.l 
JINC aO 
JHINalABSI Jl 
JMAXajMIN 
GO TO ZO 
|0 CONTINUE 
JlNCai 
jHINaMlN 
jMAXaNAX-l 

IF ( MlN.EO.O) UM|Na| 

IF ( MAX, EC. 0) JMAXaN-l 
IF (I.LT.O) ILOGa.l 
I iNCaO 

IMINat ABSI I t 
|MAXa|M|N 
ZO CONTINUE 
CONDEP'0.0 
no RS IIa|M|N,|MAK 
00 NO JJajMlN.JMAX 

IF(CONVAL.LT.AHtN|(A(tl,JJ|.A(l|*t |NC,JJ*J|NC) I 
I.OR.cnNVAL*GT.AMAKl(AMl.UJl.A(ll»i I NC t J.l« J I NC I > I 
Zgo to ro . 

IF iilog«lt.oi go to aO 

CONOEP* I ICOnVaL-A » I 1 , jJl )/I Al I !♦! INC , JJ*.||NC l-AI I I , JJ|| | 
l«FLOAT(tl*I|NC*JJaJlNr) 
return 
30 CONTINUE 

CONLOGaALOGlOlCONvALl . 

CONOEP*(<CONLOG-nLOGIn(A(If,jjll|/(DLOGln(A(||«I|NC.JJ«JtNC)l< 
lOLOGlfll A( I I , JJI ) 1 ^♦FLnAT^ I |a| iNCJjajINC » 

return 
ro continue 
rs continue 
return 

END 
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main 


»»0OC5IM fi 

Solution or oirFusiON p^oolfh rO 0 silicon on sa^pniNc 
••• NO0NAtf2£O solution ••• 


• • 
• • 
• • 
• • 
• • 


OaTa is NrAO rNOH oFC»( IN FOlLON|N6 SCOUcnCCI 


•• 


•• 


Lr iLf 

I^LOT 

ircad 


flWST C»"0! 

LIST - « Of time steps to be SKIMPfO «MlLr P»INTIn6 

iriLE - PUT 1 TO write On File *Mn also Ti> locate 

contour ROciT inN . contour FRONT MOvEMCnT DATA IS 
written on Unit 9 and concentration profile is 
written On Unit IU 

a OF time steps to rE skipped RHIiE writing on FIlE •• 
put I TO PI OT profile IN PRINT OUT •• 

put I TO read OATa from file •• 

PUT 7 to RfAD |0n implant OATA ANn TO PR 

00 redistribution. •• 

JSTeP - a of OaTa steps To rE reao from file IF,IREAD ■ I •• 
IMAXI.JHAXI - a Of 6P|n points in X and V directions aa 

RESpfCTivELT.CHCcK dimension before cHaNcInr aa 
format FtrLO . 8110 aa 

• • 
• • 
• • 
• • 


JS «"! M?'6 ,r Const, source oiff. is desiRCo* 
- PUT I IF RFoISTrIRUTION is desired. 

IF twO-STEP oiff. is desired 


orinT 


PUT> 

PUT 

PUT 

PUT 


I 


I FOR iOO cRTSTaL orientation 
7 . 

3 


FOR iIO crystal orientation 
for til crystal orientation 


AMRNT - PUT I FOR oRY 0*T6EN 
« PUT 7 FOR sTCam 
• PUT 3 FOR nITroGEN 
FIELD 3| in 

Third card: 

CSUR • substrate OiTPlNR/l .Ets 
CS - surface CONCFnTRatION/| .E i8 
Temp • temperaturf in dE&. cent* 

TmaX - NORMALiZATinN TiMe In SECOND 

this mas no effect if LaMDa is specified. as OaTa 

XolST. YDIST - width and TMICKNESSIIN MICRONI oF tmE two 


OiThIC 
oflt • 


OimeRSidNaL region |N question 
WIDTHIIN mICRdNi or The nilOE IN THE REGION 


normalized time step 
format FlrLO • BFIO.3 


imTtpe^* Specify ttpl dF impur|Tv §y putting n 
NO SPEC. |S nfCESSaRY if it is BoRON 
IMPUTT - put BOROM . ARSfN IC » PmoSPmDROUS or ant 
other Name* 

Ea » ACTIVaTiON ENfRGV of The diffusion 


on p 


01 


IF BLaNX and rOROn diff., 
OlFFUSIWITY const. OF THE 


Data Is supplied Internali v 

. , . _ Impurity . 

IF BLaNX aND bOROn oiff., DaTa is supplied internally 
FIE l 0 At,MAH.7F 10.3 


fifth card: 

ID - identification comment to be printed on top or 
profile print out 

ITE5T - PUT 0 to RfAO t*MOA FROM DATA OECx 

CSTOP - CONCENTRATfON/l .EIS aT WhICH SIMUiATION STdPS 

when The LfFT fNO corner of silicon and SaPPmIRE 
interface rEachEs This valueiouring preoepi 

FORMAT FlrLO • I 5 AN . I S • F t S • 9 

SIITh CARDIUSE IF |STp>,3i: 

ROTEMP - redistribution TEMPERaTurF 
ROTMA* “ redistribution NORM. TImf 

redistribution final TImf is i. 

HOOLT - RtOlSTRlBUTlON NORM. TIME STEP. 

XOa • REOISTrIbUTIdN initial OXIOF THICKNfSS in Cm 
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iMCwr su«r. conc. * 

*0R • ReOlStNldUTtoN InITuc OXiDr 
•MCRC SliRF. .Mas tmc Tmi 
CH 


s fS 

THICKMiSS in ch 

, - - I oxiDr 

«FNFRATro internally if 

IMPURITY IS rORON ANO NO Vai UC IS SfVEN 
fielo AFIO.3 

SrYFNTH CaRoIUSE rHEn ITEST»0|I 
LAmOA • l8m0a»»2/| .C-3 

format FIrLO • F|0«3 


fJS-lS-ElSUVflN'Si)!'*'''” 

CmaX . HAX, CONCCNTSATtON 


Rp • RANSe Or distribution* HEAN value IIn micron I 
drr • straggle* stinoaro deviation. UN Micron i 
volt - implantation enrRgt level in KEV 
field 3EIS*4,fI0.S 


• • 
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• • 
•• 
• • 
• • 
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IMPLICIT OOUBLt precision lA-MtO^Z) 

integer ORINT«AHrmT 

DOUBLE precision irB»NftLAMOA 

COMMON /TRI/ A(8R)iB|ARl»0fAH|»E(8R>«VC8Rl 

COMMON /CON/ ceilAHtAR) 

DIMENSION X(6rT« CBfAS»*€||« CB2l8i|.iRI« GCaR.ArI* |D(IsI 

DIMENSION iHTYPEftU fMPUTrIR)* IRRONIRK 
DIMENSION XoLASTiaRI. X0C4r|, DRRXO ( AS I b vP I AR I 
DIMENSION |HRTl7|tlORNT(t l•toOXrE) •INtTl2l»|STHCf I 

define computing niNCTfONS 

!!?!8!8:;t3;;4?i£7?:?l3iT5?S:s:!5'*"“'“*'''’ 

rU|RATlolROO?l*| (RAT| o*ROOtI*«I?/ROOT 
CALL ERTRaN (VilnATE.iTlHEl 

RfaO simulation INfTUi DaT* 

READ 340, L I ST • ir |LE ,l r ILC . IPLOT , iREAOfUsTEP • |M ax I , JNaX I 
READ 340, JSTP,0R|NT,»MbnT 
|F( IREAQ.E0.2i JS»E*I 

READ 370, CSU8,CS.TEMP,TMA*,X0lST,T0IST,ft«THiC*DELT 
READ 320, iMTrPE.fMPUTT.EA.OI 
READ 380, IO.ITESTiCstOP 

IF IJSTP.6T.0l REiO I Si 370, ERR* 3 1 0 1 ROTEnP ,RoTMai .ROOLT , « 0A , X«R 
IF UTEST.Ec.Oi RrAO |P,370,ERR*3lOi LAMnA 

OFFINF data: 

KP • B0LTZMANN»S COMSTaNT 

Jl. |M • a OF ALLORarLE iterations 

OvLIN • ALLOrabLE fRRor in convergence 

TinaX - NORMALIXEf) T|mf AT RHICH SIMULATION STOPS 
Ea8 and 018 - BORON OIfEUSIvITY OaTA 

data EAB,DIB /3.S9D0,i* I 700/ 

data IBRON /NHBOROi^HN ,rh ,RM / 

data NTP,|PT /IHN,IHP/ 

data kb /8.AI40-S/ 

data JLIM,0VL|M /aO, 1.0*8/ ,T|Nai /1. 00/ 
data I00X/4H0RT 0iri4HY4CN / 
data INIT/AmNITRDG.AmfN / 
data I5TM/4mSTEan / 

IF < IMPUTY.f 0 . IBRdN) |MTYPf»IPT 
IER*0 

IF I jSTPjtO. I . and. IRFaO.EO.OI lERat 

IF I IMPUTv.nE. l8RnN.AND.CA.ee. O.OOl IER *1 

IF ( INPUTY.nE, IBRON.anO.OI.EO.O.DOI ier»i 

IF (IFR.ew.ll print nsO 

IF 1 IER*ee* U 60 TO 3|9 

if (IMTYPf.FO.IPT) ITyRe»0 

IF ( IMTYPf .fO.NYP) ITyPE»I 

IF ( IMPUTV.fQ. IBRnN. anO.EA.EO.O.DQI EA*EaB 

IF I IMPUTy.F e. iBRnMtAND.DI.ee.O.nnl OlaOlB 

PRINT NIC, 10 

PRINT 330, IMPUTY.IMTyPE 
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|07 

||* 

I?? 

\n 

I** 

l*s 

I** 

1*7 

I** 

200 

in 

i§^ 

20S 
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2M 
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217 

21* 
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220 
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|Ih 

lit 
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III 
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237 

23a 
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c 
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c 

c 
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c 
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c 
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c 

c 
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NaMttlST /auT/ LIS7« ir iLCtLF : te. l»LOT,|Kf AOtjSTEa* iNAVi t.lMAM . JaTP 
I.OaiNT.AHANT.rSyB.CS.TCMPjTM** .X-iST,V0UT,0*THIc.0CLT,E*,0I»|TrST 
2.CST0PtLAHD«,R0TrMP»AnTHAX,)(0CLT,x0A«>O8.CHiJLlN,BVL|M,T iHaX 
tRITE lAtRUTI 

»CrORMAT|ON Of COHPUTmfi, fifOMSTlItC ANO RhySIcaL PaRaMeTeRS 

tf IJSTf tCOtOtOR* JrTP .«;T«1 1 ICONOal 
IFlieONDtEO.I) ANrNTbi^ 
irtAHBNT.NE.S) XOfSTaiOlSTal.O-* 

If ( ahbnt«ne.3i vntST«vD|sT«t to-a 

If ( AH rNT*NC*3) OXTHIc-OXTHIC* 1*0-4 

jMAXajRAXt-l 

KOTIaO 

INAXalMAXt-l 

|MAX2«INAXI*I 

IF ( IMAXI *gT.3) go TO R 

OELXa**DO 

If TaNRNT*NE. 3) DFI.Xa4.D><l 
. 60 TO * 

» OCLXaXOIST/rLOATl iMAX 2»3) 
k OELYaYOIST/ftOATljRAXl 
CSURaCSURal ,0|S 
CSTOfaCSTOfaltOlS 
lAMDAaLAMDA#l*D-3 
CSaCSaltOie 
TaTEMP*273* 

RINOOaiXOlST-OXTHtCUl *0-S 

GMaCM 

CRal.O 

TMXaTMAX 

If IAN8NT.E0.1) IMrTI I t«IOOXI I I 
If IAHRNT.EO.I > lMATt2)a|00X(2l 
If I AHBNT*E0*2I |MRTa|$TM 
If IAMRNT.E0.3i IMrTI I t"|Nir ( 1 1 
If |AMRNT*EQ.3) lMRTl2|a|M|T(2) 

If f ORINY .EO.t I lORNTatlOO* 


If (ORINT.ro. 21 
If (0RINT.EQ.3l 


inRNTa* I Ip* 
lORNTa* I I T* 


Calculate distancf x anO y 

00 10 |a2.|HAXl 
X( I )•( |-2 )*DELx 
| 0 IF (X( n.LC.RlMnO) ira| 

. no 20 Ja| t JMAXI 
20 V( J)a( J-l I*0ELY 
00 30 |a2,lMAXI 

30 XOLASTI I laXOA 

no 31 jaJHAXI.l.-t 

31 Vf (J)a( JMAXI -J)*( I */n.OAT( jMAxI I 

SPECIfT PREOrf CONniTioNS 

00 40 |a|,|MAX2 
DO 40 ja|, JMAXI 

CM I • JiaO.C 
CM ( I » J|aO*0 
40 CB2lt»Jia0*0 

no 50 |a2,X 

CM I .JMAXi laCS 
CBI ( I • JMAX I |aC5 
50 C82( I. JMAXI |acs 

TIMEaO.O 

1 NaO 

specify initial PRoFIlf 

IF ( IREA0.E0.2l CALt Sllg | ON ( I M A X | , jM AX I ,IC , Y , YO I ST , cSTOP I 
IF ((READ. CO. I) CtLI- XYZ I I MF I . JMF I • X . Y , l M , T I ME t K A • OF A i TmP . T I MX . J4 
ITEP.XOLAST .VOIST.aMBNtI 
IFllRCAO.E0.il TlMAXaTlMC*! .0 
ST IMEaTIME 

IF ljSTP.EQ.il IfrtNOa? 

NCAO contour FRONT MOVfMENT FROM oaTA FILF IF IREaO * I 
IF IIREAO.EQ.I) CaCL fRONT I I M i . jh I . L A t Kr • TMPR . OF b • OL T .L m . I , KnT t . x 
IF (|RE*0«E0<l) print *30. LM , T | he • TMP , OF A . T | MX . i mF I , JMF t • K A 
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288 
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307 

308 
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318 
315 
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317 
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IF (|KEA0*i0«D P»INT 81*0, TMr . ,f)r8,L6* iHl ,JH| 

IF (IReA0.NF*|) SO^To *9 

IF liMPl*NE.|NAX|.0R..|H|F|.HE*jHtit.08.KA.NC*KI lE»*l 
IF (|Ml.NE*|HAX|.0N.jMl*NC.JHAA|.08*KC*Nr«K*08*LA*Nr*lM) fE8a| 
IF l|F8.E0*|l FRiiiT i|«0 
IF I ICF.eO* I I GO TO 3(0 

start simulation STIP 

*0 IF (ICONOtCOtl) CO To Tq 
XOT IpO 
TEMPpROTEmP 
T>TEMF8273* 

DELTbROOLT 
THAXbROTNAX 

IFI AMRNT.NE.3l OFI TaorLT«TMAX 
IF( AMRNT.NE.3I T|mAXbtNaX 
^ TIMEbO. 

70 CONTINUE 
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Calculate ni 

TIMAXbTIMAX.OELT 
CALL PARAM INl.Tl 
OFbOI*OEXPUEA/(Kr*T|» 

IF IITEST.Nf.O) nrlB0F8TMAX/( (TntsT*1.0-8|B*21 
|FI AMRNT.NE.31 OFt*OF. 

IF IITEST.ea.O) OpIpLaMDA 

IF llc0N0.EO.il print 380 

IF (tC0N0.E0.21 PRINT 350 

PRINT 380. Nl.OF 

PRINT 820, IMPUTT.TEmr.oFI ,TM*X 

PbI0F|*0ELTi/I0ELT8»21 
0b(0FI80ELTi/(0Ei X**21 

PRINT initial profile 

CALL OUTPUT C X , Y , t M ak |# JM AX I , K . LM . JJ , T I Mf . TO |ST , i o » I T I ME *XOLA|tI 
• PTIMF.OTiMc.aMrntI 

IF (IREaO.EO. 2. ANn.lFlLE.EO.il CALL ARC ( I NAX I . JM aX I , | . v , L^ . T | Mr , « 
1 .OF I .TEMP.TmAX ,KOT l .X oLAST .YO lST, aMBNTI 

AMaUO 

CCb I .0 

RRat.O 

IF (ICOND.NF.ll call oXOaTa ( aMrnT .OR I NT . T .BR . 
tCC.AM.KBl , 

IFI ICONO.EQ.a.ANO.GM.rQ.O. I CM bam 
IF( tCONO.EO.a.ANO.GM.AlE.O. 1 CM BgM 
!F( IC0N0.F0.21 pRiNT 821, I MBT . | ORnT ,BR .CC .CN 

Start time step loop 

80 LMaLH^I 

STORE N TH. result for R.h.S.,R1L( NOT BE CHANGED oUR|n6 |TER. 

T |MEbTIME*DCLT 
IFl ICONO.CG.Tl STfMEarlME 
PTIMEbTMX»ST|MF 
?F( ICONO.CO.2l PT|ME«tRX 
OTIMEbRoTmaXpIT |mc*ST|NE1 
IFI AMANT.NE.31 OT»MEbtIME 
DO 80 |B| .IMAX2 
00 80 JbI , JMAX I 
80 CR2( 1 , JlacBl ( I , J1 , 

IF llC0ND.C0.il AO To IIO 

Calculate oxide thicrnpSs 

IFI AMRNT.NE.3I 60 TO p5 
DO 88 |b2,|mAXI 
88 X0( I IbXOLASTI I I 

AO TO 110 
85 PSBRfltOELT 
OSbRR/CC 

00 100 l•2.|MAX^ 

X0(11bX0LaSTI| l*PS/f ?.*X0LA5T ( I 1 PCS I 
100 0RRX0( I iBfXol I I>XO la iT( 1 I I/DElT 
XOt I |aX0(2l 
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ifo 

321 

322 

Hi 

Hi 

327 

II! 

Hi 

33 H 

Hi 

HI 

i’o 

3*l| 

3*‘2 

3<I3 

S<MI 

3 *(S 

3 H 4 

317 

3 lt 

31 » 

3&0 

SSl 

3M 

353 

351 

35S 

!H 

3*0 

341 

342 

343 
351 
34S 

344 

347 

348 
341 

Hi 

111 

HI 

378 

ill 

381 

382 

383 

38 ^ 

38S 

384 

3«7 

388 

388 

3?I 

311 

314 

317 

318 
311 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 

c 


c 

c 

c 


ir 1 1 v6IST-o.is**(h I n .tT.o, 1 p- i ) r,o to 3io 
llO rONTiNUe 

OCLvall •/7L0ATUMtX) )•( YOIST-OtHStXOLASTiai > 
tP(AM8NT*Ni.3) Pa ( OF | •DEUT ) / ( 0 ElT**2 ) 

JJ«Q 

jHaO 

STA8T iteration LOoR 


120 jjajj«| 

IF IJJ*g7*«II|H) paint IQOi CTeR 
IP IJJtGT.JLiM) (in To 2*0 


Transfer soln* vector tRom last iteRa* 


for cal* of 6 


00 130 I>1*|HaX2 . 

. 00 |30 JaltJPAIl . 

130 cRi I •JiaCRi ( t ,«n 


CaL*6 

DO IRO I>Ii|HAX2 . 

DO |10 J«I*JHAX| 

RA>RAT|0(C8I I ,J> ,r5U8,Nl 1 
R0«R00T|RA) 

110 Gl 1 »JI«FU|Ra,rOI 

SOLVE STSTEH in r oIREfTlON 
DO 220 1>2 *|MaX| 

Calculate the coeff* in y direction 


151 


ISO 


SR*0* 

00 ISO J«2*JMAX 

G1>IGI I , JI^GI 1*1 ••ij 1/2.0 
62>IG( I ,J)*GI I »J*| 11/2,0 
6S»(GI I . J)*GI 1-1 , jl 1/2.0 
GRXGI I.JI*GI|tJ>l>)/2.0 
|F(AMRNT.CQ.31 GO to |S| 

SR*0.15*0RRxOI 1 UfTPljl-l .01 *OELT/OELY 
CONTINUE 

At J)«*P«GR 

R(U)*1 **a*(GI«63|aP«IG2*6N>*SR 
0( JI"-P*62*SR 

EIJ1>RHSI0*6I .63,rB( 1-1 , JI ,CA2( I t J) tCBI UI *J| I 

Put rounoanv condition On t axis 

M?|aAl2UBI?) 

AI2|a0.0 

IF (1C0N0.NC*11 60 tO l7o 
IF ll.GT.KI GO tO IaO 
El JMAX|aClJMAX|.OIJMXXI«CS 
0|JMAX)*0.0 
GO TO 180 
CONTINUE 

81 JNAXlaB UMAX 1*01 JmAXI 
OIJMAXlaO.O 

IF |AM6N?,E0.3| ORRxOiIlaO. 

Ha«-DEL V* ll./CM>0.<«G> *ORRXO ( 1 ) 
hR«2.»0FI 

HK |aHA/ IHR*6l 1 , |RAX I 1 > 

HK2*MA/|HB*GI I .jPAX I 1 

B(jMAX)aBUMAXt*<<l.*HK2l/(l.*Hrl))*0l|MAX| 

0( JHAX laQ. 

180 CONTINUE 

CALL TRIOaG I2,.|MAX» 

convert matrix solution 

00 |10 J«2.jHAX 
CRI 1 1 . J»>VU) 

113 CONTINUE 

PUT roundart vaLUfr In t 
CRI I I , 1 )"CR1 < I ,7> 


140 


170 


< 
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<tOO 

<«0t 

<|02 

•tOI 

•tO* 

<i07 

«OI 

’f? 

HiZ 

HlS 

<«u 

mi 

MIS 

MIS 

MZO 

M2| 

:i! 

HZi 

427 

42t 

424 

420 

\\\ 

424 

:3i 

ill 

s3f 

mso 

MH 3 

MM 5 

MM 4 

ms; 

4%^l 

4i0 

4S| 

: I 

mSm 

mss 

mSS 

m57 

MSS 

MSS 

MSO 

MSI 

MS 2 

mS3 

MSM 

MSS 

MSS 

MS 7 

MSS 

MSM 

M 72 

;» 

J?2 

477 

478 
474 


200 

2l0 

220 


ir ( icono.nc* t > 60 TO 2in 
ir f 1 . 67.1(1 60 ro 2<|0 

CAM I .JHAll )>C4 
60 TO 220 

COM > .JHAxn«C8l M,.tMAX) 
60 TO 220 
CONTiNUe 


COM I .JHA*M»M« •♦Hk 2»/M .•Hljl M4C0M I . 
irTcoi 1 1 1 JNaXM *14* t .n*200i cam i • JHAXI )•! *0*200 
CONTINUE 

Put roundaov vaLueA in axis xiACotooic oounoaovi 

00 220 jal ..IHXXI 

COM I .Ji'COI (2..|) , 

220 COM IMAX2.J)>C0| UMtX.j) 

check for CONVCPeCNCC 

ir (JJ.CO.M 60 TA 12A 
|CK»0 

00 240 !•! .IMXXI 
00 240 J«I.JHAX . 

ir (CBMI.JI.I E.O.Oi 60 TO 240 
CTEST»0*BSMC4l( I.J»-CA( I .JM/CBM I ,.lM 
ir ICTeST.LE.AVUlMl 60 TO 240 
CTER«OH*X|lCTrST,rTC(n 
. 1CK*I 

240 CONTINUE 

IF MCKtNl.O) 60 tO 120 
jMb JH^I 

00 241 laltlNAXI 


ORIGINAL PAUL lb 
OF POOR QUALITY 


00 241 jBl.jNxXi 
241 IFICAM t*J).I.T.O.AO) 
IFIJH.CO.tl 60 TO I20 

PKINT RESULTS 

|5bLM/L1ST*UST-l M . 

IF MS.E6.0) 60 To 2SA 


cot ( I .UIbCBM I • J| 


2S0 


IF (TINe«6|*TIMAX,ANo.CS.NE.0.0) 60 TO 2S0 
IF(CBM2.|T.6C*CST0P.AN0.lc0N0.Ca.l) 60 TO 200 
IF(CBM2tM.6E*CST0P.AN0.CStC0.Q.A.AN0*U4TP*e0*M 60 TO 2&0 
GO TO 2A0 

CALL OUTPUT < X • V , | NAX I » JN AX I . K i LM . J J . T I Mr * T0| ST . | o * I T I MF . XO , 

. JMaXI ,V,T|ME.OELT»O tLT.TEMP.TMAX,ir ILF.fTY 

IRE. AMRNT.X0I2M 

STORE TRANSIENT OatA 

2AC IF MFlLE.EO.OI 60 To 270 
|T>LM/LF|LC*LFlLc*k^ 

IFI IT.E6.0I 60 TO 2A$ 


IF(TINE*6E»TIMAX» 60 Tp 2A% 

IFtCBI <2. 1 I .6 E*CStOP.aND. I rONO.EQ. I I 60 TO 2AS 


:»:*«»? i»t«T 


locate CONTOUR POSfTlON AND STORE OATa 

2 A 6 continue , ... 

CALL front MMAXI .JMaiI ,lM,K. temp. OFI , DEI T.JSTEP.OiKOTI .XO lSTi 
270 CONTINUE 

60 TO next Time stfP 

IF ICRMT.M.GE.CSTOP.ANO.ICOND.FO.M go to 240 
IF ITfMc.GE.TiMAXt 60 TO 2TD 
XOI M•X0I2> 

no 280 l■2.tMAXl 

280 gX^^AgjM >•XOM I 

240 IF ( JSTP.gT* I • and. I COnD.eO. I I GO TO 300 
. CO TO 310 
300 T|MAX«T|ME*I. 
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•Sim 

:!! 

Mil 

MM 

M»0 

Mil 

:u 

SSI2 

mm* 

100 

lOl 

i8i 

•8> 

10 * 

107 

SOI 

SOM 

II? 

%tz 

il3 

IM 

SIS 

SIS 

It: 


IF I *MRNT*Ne.3l T |M*X«T IMC^MOTh* . 

IC0N0*2 
. GO TO *0 

310 MSo 

iiS JIShaI ltM6!loii?*}HfuilTY - »,**»M,2*,MM|iUIITV TVMt • *,»Ml 

|M 0 roi«»T ( iNe*208i«*«* mg^ocmosition Cycle •••*) 

350 rOIMiT (IhT.20*,«*** mEOISTHBUTiON CYCtr *•••} 

3*0 rOIMXT (IllOi 

ill roiM»T j)/HMOifoo*NI • ;o>0t3.s«i«ori . •lOlO.s//!^ BIO Si 

I ON SAl*MH|RE •••*//|h0»3|I»*N0Mm«LIZED MOLUT I ON* // |m *IOX*I&*M//i 

1 . IMTMINSIC OlfrUMlMlT^ OF » • I * * M* M i IS t * *T * * «* t f * 0» 3 » 2ll • • oE4» cEiiT« 

2t//iH .l0Kt»L*«0*-" L*M0**»2 • 0FI»TM*K/tYH*K«|.C«M**2l*//lH •Ink. 
3»L*ND* • t,nl0.3,S«,»r0R N0IH*LI7*TI0N T»M| ■ ••Cl2,*,2X.*SCC.*l 

•IHo’si ••it6M*6*^?nN*C0EFF» • 'iCiioO) 

M30 format IlHO.IOXt**** initial MlOFlLC aT tine step • '•!Gi2||»*TImE 

• tOIO*lt2X,«|MAXl • ••II0«2 Ii*JM 

. ARE obtained fGOH ThF OaTa . 

I . file M .inX.tTtMP ■ '.FIOtX/lH «ieX»ttAM0A«*7 • »»OfO. 

23/IH .lOX.tTIME MYER • ♦,MO/|H .IOX.MMaXI ■ ‘.IS/lM •lOXt'JMAtl 
3 • ^ * f S ) 

MSO format I ImO,IOX,»***OaTa InRUT FRROR ♦••♦/lMO**X,«RON$TRf am TrRHiN 
1 ATEO* I 

END 


3 • ,0I0i3*2X/IH .IaX.InORMi TIME ■ • 
Maxi a *, 110. 2Xi«oxlOE GRID a ••!*) 
MMO format UMO.IOX,**aa fOtLOwINS DATA 


BPRT SOS7.SmIEJNi.MAIN*RL0T 
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>t|CR0N«SoS2( i I tSHNr JN 

sumiouTiNt sHnejtuCSiiiiiTDtsT .jmaxi .x.TiHr •oiiTtOei.v. 

• TC»1PiTM«X. IFlLC. ITYpC, 1 AMBNTiIUI 


s 

7 

10 

U 

13 

u 

14 
17 

II 

li 

ii 

21 

2 * 

27 

21 

I? 

II 

37 

31 

39 

90 

91 

92 

93 
9H 
9S 

94 
97 
99 
99 
SO 

li 

15 


SI 

S9 

{? 

42 

43 
4M 
4S 

44 

tl 

*9 

70 


implicit pOUBLC l>«CCI9l0N(«*H»0»2l 

COMMON /CON/C(*N»At) . 

01 mens I ON CNCTUMI iCMI i*«0 ,XI I ) «v(*9| 

OCLavOiST/FLOXTl jmAXi.I ) 

tr I IaM 0NT«NC»2) OFtaort/VDtST 

V*J«0*0 

VNjIaOtO 

DO A laliJMAXI 

I JaJMAX|-r*l 

VIMaVDIST-XltJ) 

IF 1 1 «MBNT*NC*3I Vllla(l*ll*l|«/rLOATIjMAil«in 

jNlal 

ICKaO 

00 6 jajMAXi 1 1 
|ajM*Xl*l«J 

ir iclitiitfTicsuii .AN0.ICK.CQ«0» .iNl*l 
|F(CBim*6T*C$UBI ICK*> 

IFICBI ( n*«T.CSUM) JNal 

IFI JnI ! «(JMiX|*l|*Oll.jN.e0*jMAX|) 60 TO 7 

YNJalvOlST/l JMAll.l ) I • < CO I I JN I •CSllB ) / I CO | I JN ) aCO | ( JN« | I ) 
IFI t*M0NT«NF*3> YHJayiKJ/fDlST 

r«J|aOEt*ICBUOaCMl Oni ) i/ICBl (JN|*I l-CBl l«INI 1 1 
VaO 

SifiMAaOtOOO 

0a0*00 


STaMT INTCGMaTiON 

IFlK.fiEf JMAxn GO TO 9 ^ 

OaO^ICBI IK«| )«CB| (K) |•0EL/2•0 
IFlK.LCtJNl .ORtK.ftEt JN> GO TO I 
SIGIaGIOABSlCNETlKi I . I TYPE I 
9|G2afilOABSlCNETl««l I t ilTYPCi 
SIGMAaSIGMAal YlK*t laYlKn*! tSIG2*BIGI 1/2.0001 
GO TO 1 

continue 

IF(JN.E0*JMaXI ) VaJao.O 

IFI JNI .EO. I .ANO.rilM JN> ) tGT.CSUBI Y«J|aO.O 
BiGAaGlOABSiCNETl.iNai t I • ITyME ) 
BIGBaGIOABSiCNETi |N| ) | i ITVME) 

SIGCaGlOABSlCNFTl |NUt 1 I , ITYME) 

VJUNCaTI JNI^YBJ 



VWJIaVIjNI »♦YllrJ| 

^IGMAaSIGMAA'SlGIVj 
BSal .eOO/SlGMA 
NII|TF|At203| 


: • • 


"S,G. 


„.MTF|A.203| Tb j I I Y ji(NC • I V a I V. 

20) fOBMATI/.&X, •junction is at <«2niG*9.« CM*tBX./t 
•GX.*ShEET MfSISTamCE a •iOlB>9»3l« 

• /3X.*INTrGMATE0 I MPu* I T Ya • , 0 1 5, 9 1 


201 


»BtTE|A.20l I JNI.jN 
FOPMATI *0« iSX,* jNa*i2tB/| 


AF ACa| .0 

IFI lAMBNT.NFtSl Ar AC a ( YD I ST -0 • 9Ga 1 0 I / 1 • 0*9 
V JUNCaYjUNCaAF AC 
•SaBS/l AF AC* I .0-N I 
OaO*AFAC* I .0-9 

IFI IFILE.FG. I I CAi I SmF HE ( T |HE .O ltT.OEtY . tempi 
• TMaX.BS.YnJI •V.iiiNC ,n I JMAX I ,VD|9T|X6.IAmBNT| 
RETURN 
ENO 
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MICR0N*SnS7l I I •mA|«|-^LOT 

•• iSOrONCf NTH»T ion PloT PP 06-»*H •• 


I 

! 

n 

I 

4 

7 

• 

f 

10 

11 

12 

u 

IS 

u 

12 

17 

it 

23 

2 «i 

II 

it 

i! 

m 

3S 

!t 

u 

if 

fo 

fl 

72 

7S 

2 t 

ft 

7f 

11 

ll 

IS 

§7 

ftt 

Sf 

40 

4| 

t! 

4<l 

44 

44 

47 

4t 

4f 

;? 

72 

» 

n 

ii 

7f 


3 

C% 


ff 


40 


C* 

c* 


c* 


r|t«7 CA 




IN r»OP OAT* OCCK 
U - • Of Tint STfP< TO tc 

*l-Af»„lN,r80M 0*TA rilE, 
SECOND CAKOI IIT * 


- T|«E STEPS^TO IVE SK1Pp2o^»huC 
PtfiTTlN# PSOffLE* DEfAULT PLOT AT 
T|HE STEPS I AND *T iL. fT 


PIELO 13* 


01 HENS I ON CaU<«tAA> tlrUNIiTCiAA) .AO(AN) 
OATA IMAI.ThAI/S.O***. 0/01/3*4/ 
data tOOLiiSLNK/iHSiiM / 

OATA M6T/0*Ot7S/ 

PORMAT 113) 

• ^RCAO OATA PILE I I 
00 Too La) ILL 
PCA0(7|N) NaRK 
PORMATI INOiAA) 

IPINARK.NEiIOOL) rO To S 
REA0(7.n DPI iTmav iTemPiYOIST 
REA0(7i2) K.IHAXi.JNAllitAMSNT 
IKalMAXI 
jKaJHAXI 
PORHATUH 
PORMATI IH 


ORIGINAL PAGL 1 ." 

OF POOR QUAUTY 


•SElSitl 
iSI lO) 

RCA0(7il) (xCm,|a2||K) 




REAOlTill lYCin,ta||jK) 

REAOI 7 iN) Mark 
rONTINUE 
REAOI 7 i 2 l . 

REA0I7|I) T|Me 
00 ff JaJKi I 1*1 
REA0l7tll icon i.n ilapilKI 
00 40 fat.lK 

IPICnilf JI*LCtl*P-30l CBI l,J)a|.r-30 
CONTINUE 

RCAOlTii) IIO(n,ta 2 .?K) 

ISat/ltT*| jT»L 


|PIIS.E0*0) 60 TO II 
IPIL.EO* l•OR•U•re.LLl 60 To 
. , «0 TO 100 
II CONTINUE 

• INITIALIZE PLOT. 

CALL PLOTSI iQ.OilOiO) 

• RESET PIN TO 0R16IN 


I I 


5 ‘il. 


C* 


0R16IN 
PLOTII ,0*Z,0.-3» 

• ORaR roroer 

CALL PLOTI0*OiVMaii2i*1) 

CALL PL0TIXMAX|VMaAi2,-| ) 

CALL PLOTIXmaXiO.QiZ..} I 
CALL PLOTI0*0i0»n.2i*| ) 

• ORAN OXIDE MASK 

XPaJXC ( IKI-P. I . , . 

IPI i AMBNT«Ne*3l Xaa| Xc I I K I -2 * E *4 I / I • E-f 
OXa|R,/XPi*XCIKt .. 

CALL PLOT lOx ,yM** , 3) . 

call plot jOX iYM»**0*2,2i-I > , 

CALL PL0TiXMAXiYMtX«0.2tZ*>| t 

call plotiimaxi ymax,2,-h 

xaXMAX 

PATat,-OX 

JNCalMT < PaT/ 0*2 I 
00 10 lalilNC 
Xax*0*Z . . 

CALL PL0TIX.YMAX*O*2,2i-| I 
CALL PLOTIX.YMAX.S) 
lO continue 

call PLOTIO.OiO*n.3l 

• horizontal 6R|n LlNfS 
no 20 l»l.f 

Yap.N.PLOAT ( I I 
I NNTa*3 

IP 11*16.4) LNPTa*l 
CALL PLOTIC.O.Y,!! 

. CALL PL0TIXMAX|Y.7iLN«T| 

20 continue 

• vertical 6RI0 It NEK 
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00 )C taltlf 
laOt'l'rLOATl I ) 

1 NItf . , 

If (MeD(l.il*rO*Oi kN«T«»| 

cSut 

•®Noi|HoNT4t *«M NU«*C«1 

n|55"”- 

iflO|l*LT*e.i I lc*ll 
ir(o|i*kT*o.h is>f ... 
ir ioli*kT*o*oi ) nKX«i»/ie* 
in UMiNT.cotsi 60 To Si 
| 6«l 

)t rONTiNue 


kaO*0«rLOATl l>{ 

y«»o.a 


CAUU NUOICRII»T»M6T*V«WiO*OiU> 

"** t^veilTtCALtAiis niiHie»s 

in l«M|NT*Nr *)> V«KaVc<«IK»/l«r*A 


OCVaTAK/10* 

00 to . 

VAlaOCTartOATl l-l I 
|aa).0*H6T*0*l . . 

VaVMAI-0*0*ri0ATI l?» UO 
CALC NUMtCNi I t V |M6T| V«t 
tC cONTINUC 

rALt •|VHtOt(*)*|it.JiO.! 
rALL SVMt0LI*lta«».|«0. 
CAfC lvN»0tl<«>2tt.3*0.| 
rALL STMtOU]!*2*t.l.O.} 
(ALL tTHtOL<;*ltA.TtO« 


TAIL STHKOLiyt^.A.TtO.I ItOCiA 

IKISi;!;:8:{:|:8: :f:S:8::il 

CALL STM00L(y*0it*2i0.i i2i0«0«-l I 
call tVMtOLI*l*Ott.l (0. ! •3i9*0i*l I 
(ALL fYHtOL(A*0iA.Tt0*l •A*0«0.-n 
(ALL «VMtOL(*l*Oia.T,O.I itiOiOi*! I 
CALL tyotOL«OtO|t.'».0.2{» 

CALL tVMtOLlO* if *2.**l lO. > 

8Stt :i!Hk!i::;t;j!4:j;5i::8;:J! 

(ALL fTOOOLlOiOif.l le.li’TeHACtAT 
CALL »iUOit6(2*Ait*l i0.2|TeMPiei0. 
(ALL SyMtOLfO*OiA.tiO«2i*T|MC STr 
Tt*<«L»« . . 

CALL MUHeCtl2»Ai«t.fi0.2iTLHi0«0«> 
(ALL SyH60L<0.0|A,ftC.2i»T|Ot 
CALL NUH0lKl2*AiA.ti0.2iTlMCi0<0« 
00 200 lla|,4 

r0NVAL*l0*0««FL0ATl20at|*li 
CALL ALTCOMICtiCOaVAL.IIilKiJKl 
cONfiNue . 

(ALL alot io.OiO*o.yy«i 

rONTtxuc 

ATOA 

(NO 


t«H6T 

o*oi n 


1 1 0C2n* lOiOitl 
I loim lOiOift 
i lOCl A* lOiOifl 
I.o|iy«i5i0ifi 
liOClA*iO.Oit| 
liOClA* lOiOtOl 


2iOiO«*l I 
3i9i0i*l I 
AiOiO.*l » 
ti9iOi*ll 

• a • tOiOi I A I 

• 0* I * I i 
•|•lO•ltl 
OftiO*|A| 

•TCMAttATUAC a •,0*Cil**l 
TCMPiCiO.ni 

•T|MC STrA a *.0*0il‘*t 

TLH|0«0«a| 1 . « « .«» 

• T|Ot a •,0*0il*ll 

TlMCi0<0«2i 

*11 

•iKiJKl 



